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A. BACKGROUND 


Since the introduction of the Sampling theorem and tne 
matched  fiiter, digital communications techniques have 
developed into a highly proficient discipline. The marriage 


of this discipline with the rapidly expanding space program 


has resulted in communication satellites employing а 
multitude of digital signal modulation  tecnniques. А 
nodulation technique is a method of transmitting the 


information contained in a message by varying cr modulating 
the characteristics of a carrier waveform. These metnods 
offer a range of advantages and disadvantages depending on 
the specific characteristics of the modulation technique 
employed. The applications of the various signal modulation 
tecanigues likewise vary. In order to understana this 
exciting new field it is necessary to look at some of the 
aspects of satellite communications systems in general. 
Then an iLvestigation will Бе made into the general 
attributes of digital communications and their relationship 
to the satellite system. | | 

Ihe ability to understand these digital Signal 
modulation techniques is the first step in being able to 
intercept, identify and demodulate unknown digital signais. 
These digital signals, transmitted from unknown sources, are 
believed to display frequency characteristics peculiar to 
the mcdulatioun technique employed in tne encoding process. 
The digital computer offers unique opportunities den 
simulating these modulation techniques, in siynal processing 


and in decoding of an intercepted signal. 
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B. SPECIFIC GOALS 


It is the specific goal of this thesis to investigate 
the open literature on the topic of digital signal 
modulation techniques in satellite communications. This 
includes a basic understanding of the communications 
satellite system as well as the specific techniques employed 
in signal modulation. Additionally, опсе а basic 
understanding of these digital signal modulation techniques 
is achieved, computer code in the FORTRAN programming 


language will be developed which simulates these modulated 


signals. The statistics of the time-varying Fast Fourier 
Transforms (FFT) of these Simulated signals will Бе 
investigated. Tue purpose o? this analysis will be to iead 


to fcllow-on research in the area of signal analysis of 
intercepted digital signals whereby they can be classified 
by their FFT as employing a specific digital signal 
modulation technique. 


C. SCOPE OF THE PROJECT 


Althougn an indepth analysis of alli the digital Signal’ 
modulation techniques which will be introduced involves 
considerable higher ievel mathematics, it is not-within the 
Scopé of this project to delve heavily into the mathematics. 
It would be advantageous, however, for the reader to have 
had integral and differential calculus and an introduction 
to statistics. Also a course in electrical engineering may 
prove helpful but is not essential since in actuality it is 
the electrical engineering aspects of satellite 
communications that this tutorial is attempting to present. 
The computer programming will be accomplished utilizing the 
techniques of software engineering and top down modular 
design. Existing blocks of code will be used as nodules 


whenever appropriate and available. Again 1% 15 the 


ultimate purpose of this project to examine a variety of 
digital signal modulation techniques and develop computer 
code that simulates common Signals used in satellite 
communications that have been produced by one of the many 


digital signal modulation technigues to be investigated. 
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II. UNDERSTANDING SATELLITE COMMUNICATION 


А. HISTORY AND APPLICATIONS 


Forty years ago, in 1945, Arthur Clarke first envisioned 
the use of space stations placed in geosynchronous orbit for 
communicating to different points on the earth [Ref. 1]. 
Nine years later in 1954, Ј.Е. Pierce of Bell Lakoratories 
performed a system analysis on such a communications system 
(већ. 2]. In 1957, the launch of Sputnik demonstrated the 
feasibility of using a satellite for just such an 
application. However, by 1961 the only satellite 
communications technologies which had been demonstrated were 
the Ccurier 1B satellite, a short-life active retransmission 
teletype conmunications satellite ina 1000 km orbit, and 
the Echo I balloon іп а 1600 km orbit which demonstrated 
pasSive reflection of powerful microwave signals from one 
earth station to another. Active microwave communications 
demonstrated in Projects Telstar and Relay were years away 


(Ref. 3]. The first geostationary orbit was achieved by 
Project Syncom in 1963 [Ref. 4}. 

In 1964, communication organizations from several 
countries joined together to form the international 


organization of INTELSAT (International Telecommunications 
Satellite Organization). INTELSAT's purpose was to develop 
a satellite network which would provide truly glotal 
communications capabilities. This resulted in the launch in 
1365 oit the world's first commercial communications 
satellite, INTELSAT I, also known as "Early Bird", With 
"Early Bird", telecommunications utilizing satellite relay 
were established between the United States and Europe. 
[Ref. 3] 
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The reliability of these satellite systems has improved 
dramatically since INTELSAT I in 1965. Reliability of 
individual links in the system approach 99.99 percent or 
higher. Total system reliability exceeding 99.9 percent is 
common [Ref. 5], making satellite communications more 
reliable than many other modes of communications. Ig tkis 
sense, reliability is a measure or the probability that no 
failure will occur in a respective channel or in the system 
during the design life of the.satellite. 

In crder for the operability, capability and reliability 
of these systems to have developed at this pace, 1t was 
necessary for the technologies associated with them to 
develcp as rapidly or even more rapidly than the systems 
themselves. Engineering sciences and specifically the 
fields of aerospace and electrical engineering historically 
have required between 7 and 10 years to take a concept from 
operational reguirement to full scale operation. This was 
not the case with the concept of communications satellites 
which nas seen four generations ot INTELSAT  sateiiites 
within a decade's tine. [Ref. 3] 

The "Early Bird” satellite weighed approximately 38 
Kilograms and possessed limited power and bandwidth capacity 
ehabling С to Carry only 240 two-way telephone 
conversations. Today's communications satellites represent 
order-of-magnitude improvements in many important operating 
parameters such as fower and bandwidth. Additionally, 
increased effective radiated power from the techniques of 
stapilized earth pointing antennas have greatly increased 
the capacity of later generation communications satellites. 
INTELSAT V; the current generation ot communications 
satellites, is а three-axis stabilized platform using not 
only the 6/4 GHz frequency band (6 GHz гесеіуе/і GHz 
transmit) as in earlier generation satellites but also a 500 
MHZ candwidth available in the 14/11 GHz band. The 
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separation in the receive and transmit frequencies 15 
necessary to prevent interference during simultaneous 
operation ot the receiver and transmitter. Ancther factor 
increasing the capacity of present generation satellites is 
the increase in primary power available in the satellite. 
(већ. 31 

The technologies which have contributed to the evolution 
of communications satellites come from two primary sources, 
namely technology from the Space program of the 1960's 
supported by NASA and DoD and communications technology due 
largely to commercial and private sector contributions. 
(већ. 3] 

Electronic devices and components have contributed 
Significantly to the rapid development of satellite systems. 
These electronic devices range from something which is now 
considered basic, i.e., the transistor, to devices such as 
Traveling Wave Tube amplifiers, lightweight antennas and 
antenna feed systems. There have been major improvements in 
Satellite power sources including more efficient solar ceils 
and high storage Capacity/lightweight batteries. 
Additionally, two significant aevelopments in electrical 
engineering have made modern day digital signal processing a 
reality. They are the sampling theorem and the matched 
filter. (Ref. 3] | 

Communications satellites of the future are likely to 
utilize onboard signal processing. Sigral processing 
functions such аз signal reshaping, Switching and/or 
multiplexing and compression could soon take place on the 
satellite due to the reducticn in size and weight of the 
necessary hardware. Operation at over 100 megabits per 
second are envisioned. Onboard signal processing will 
greatly reduce the exrense and complexity of earth stations 
thereby making services of a satellite available to a wider 


range of small and gecgraphically disperse users. [Ref. 3] 
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Satellite communications hold great voromise to provide 
service to a multitude of users over a wide area. 
Applications lie not solely in retransmission but also in 
data collection from that same large geographic area. 
Military users have definite applications in these areas 
when warning, intelligence and surveillance systems provide 
digital inputs to a central source. Although most of the 
present applications for satellite communications are still 
provided by the government, commercial applications nave 
Seen tremendous increases in the last several years. 
[Ref. 3] 


B.  ORBITS AND LIMITATIONS 


Satellite orbits are Generally categorized as either 
eguatcrial (0 degrees inclination), polar (90 degrees 
inclination) or inclined at Some angle other than 0 ог 90 
degrees relative to the spin axis of the earth as 
illustrated in Figure 2.1 (Ref. 6]. Each satellite orbit 
has a characteristic velocity which is dependent on the 
height of the orbit and the orbit's eccentricity.. 
Eccentricity 1S a measure of the degree to which the orbit 
approximates a circle. A circle has eccentricity equal to 
zero since the focus is located at the center. See equation 
2.1 and Figure 2.2. 

À communications satellite in elliptical orbit about tne 
earth obeys Kepler's second law of planetary motion. That 
is, a satellite's constant angular momentum about the earth 
means that its areal velocity also remains constant. See 
колге 2. 3. Of particular interest is the satellite 
velocity at apogee (Va) and perigee (Vp) given by equations 
22221110 2-3. 

Note that for a circle, e = 0 and Ra = Rp. Therefore a 
satellite in circular orbit about the earth has an orbital 


velocity aS given in eguation 2.4. 
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Figure 2.1 Satellite Orbits 


Most communications satellites are placed ina circular 
equatorial orbit where the desire is to stabilize them over 
a fixed point on the surface of the earth called the 
subsatellite point. This type of orbit is called 
geostationary or stationary. Àny satellite which has an 
orbital period equal to the period of rotation of the earth 
is called synchronous or geosynchronous. These terms are 
used almost interchangeably in most literature. As 
mentioned for a geosynchronous communications satellite, the 
orbital period of the satellite, T, must be equal to the 
period of rotation of the earth. Kepler's third law of 
planetary motion can be rewritten to yield equation 2.5. 

The period of revolution of the earth for a sidereal day 
is 23 hr 56 min 4 sec. For that given period there 1s only 


one satellite altitude as expressed by Kepler's third law. 
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O< e< 1.0 
Circle. has 
e = O. 
Figure 2.2 Eccentricity 


eccentricity (e) = c/a (eqn 2.1) 


с -= distance from focus 
= semi-major axis 

b = semi-minor axis 

Ra = radius cf: apogee 

Rp = radius cf perigee 


By rearranging terms that altitude is given in equation 2.6. 


Since the orbit is circular, а = Ra = Kp and the height of 
the orbit above the surface of the earth is û = a - Re; ог 
35,804 kao. 


One disadvantage of a geosynchronous communications 
Satellite is the lack of global coverage. These satellites 


provide excelient coverage of the subsatellite point and 
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Figure 2.3 Kepler's Second Law of Planetary Motion 


(k/Ra(1-e))**.5 (egn 2.2) 


_(КИВр (1+е) ) **.5 | (egn 2.3) 
GMe = 3.98866 x 10**-11 п3/52 
universal gravitational constant 
6.67 x 10**-11 N nm?/kg? 
newtons - m kg/s? 
mass of the earth = 5.98 x 10**24 kg 


(k/RC) **,5 | (eqn 2.4) 


Re + h 
radius of the earth = 6.37 x 10**6 m 
satellite orbital altitude 


2 pi a**1.5 / k**.5 (едп 2.5) 


period of rotation 


semi-major axis 


Va 
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а - (T/ 2 pi)**2/3 (k)**1/3 (eqn 2.6) 
a = 42,173 kn 
laterally to a latitude of atout £809 [Ref. 6]. This is 


generally no problem for commercial applications since the 
polar regions do not require a Significant degree of access. 
The military, on the other hand, does have an interest in 
communications in the polar region and therefore has several 
communications satellites that have orbits inclined at 
various angles. This type of orbit generally has 
disadvantages of lack of continuous coverage and a much nore 
complicated system of ground tracking and receiving 


staticns. 


C. FREQUENCY BAND CONSIDERATIONS 


Although there appears to be an infinite  nuurer of 
frequencies available for communications, restrictions do 
exist as to those which are practicable. Limitations on 
available frequency bands for satellite communications are 
due tc the need to select segmerts of the electromagnetic 
spectrum which reduce noise and interference and are nost 
favorable in terms of power efficiency and propagation 
distortions. Trade offs must be made to arrive at the 
optimum frequency for a particular application since single 
frequencies seldom coffer the best performance for ali 
variables. The problems arise when consideration is given 
to the number of users requiring the same frequency bands 
including terrestrial communications networks. Since the 
problem of interference is global, a worldwide organization 
has been established to assign trequency bands for various 
applications. This organization 1S called WARC, World 


Administrative Radio Conference [Ref. 6]. Illustrations of 
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the pcrtions of the electromagnetic spectrum in question and 
current allocations of satellite frequency bands are shown 
in Figures 2.4 and 2.5 [Bef. 6]. 
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Figure 2.4 Electromagnetic Spectrun 


In general, as long as the hardware and technology will 
Support it, higher frequencies are more desirabie and more 
in demand because they offer higher theoretical Capacity. 
This 15 due to the fact that only a percentage of the 


Carrier frequency is capable or actually transmitting a 


Signal of a given bandwidth. Higher frequencies wouid also 
experience less interference with existing land and 
satellite systems. A further discussion of the bandwidths 


available as a function of frequency will be included at a 
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Allocated satellite frequency bands, WARC 1979 (f = uplink, 
| = downlink, — = crosslink). 


Figure 2.5 Satellite Frequencies 
laten pornt —under the’ topic ог actual digital signal 
techniques. 

As an additional note, in order to assist in the ability 
to handle more signals in an obviousiy limited 
electromagnetic spectrum, several muitiple access scheres 
have keen developed. The mest common tecnnigues involve 
rrequency-division multiple access (FDMA) where the 
allocated satellite frequency band is divided among the 
users into specific uplink and downlink frequencies. Next 
15 time-division multifle access (TDMA) where the satellite 


frequency band is shared by all users by carefully dividing 
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the time any one user has access. Finally there 15 


code- division multiple access (CDMA) which involves 
modulating a specific address code which has been 
Superimposed on the signal directly onto the carrier. In 


this manner all users share the satellite frequency band and 
only those receiving stations that can demodulate the 


address code can receive the specific signal. (Ref. 6] 


D. WHY DIGITAL MODULATION 


In general, digital signal modulation techniques are 
superior to analog signal modulation technigues used in 
satellite communications for the following reasons: 

1. Compatibility with digital computers 

2. Economic advantages 

3. High aegree of flexibility 

4. Less susceptible to interference 

5. Quality of signal independent of transmission 
distance and network makeup 


1. Compatibility with Digitaleeonpucers 





Clearly one of the most distinctive advantages of 
digital technigues over analog techniques involves computer 
applications. Properly formatted digital Signals can be 
used to represent any analog signal (more on this under the 
discussion of the sampling theorem). Once in the digital 
format these Signals can be easily manipulated within the 
digital computer. Arithmetic operations can be applied as 
well as logical operations. The signal can te stored 
Without alteration or delayed and theretore can be used to 
"simulate" real physical situations. 

The hardware associated with digital circuits is 
free from drift or aging and does not require calitration. 
Additionally digital circuitry is compatible with present 


day integrated circuit technology allowing a standardized 
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building block construction approach. The operating 
Characteristics of these systems employing a digital 
computer can be changed by altering the software rather than 
the hardware as is the case fcr analog systems. Finally, 
the use of digital computer technology allows time 


multiplexing. 


2. Flexibility and Economy 





The flexibility and economy of digital satellite 
communications comes from the fact that more and more 
processing can be done onboard. This allows the uplink and 
downlink to be compietely separated. This regenerative 
nature makes it possible for low error rates and nigh 
reliability through the use of digital techniques not 
available to analog systems. Because digital signal 


processing or multiplexing is less costly than for analog 


Signals, Simpler and cheaper interfaces between earth 
staticns and terrestrial communications networks are 
possitle. Additionally, there are reduced production costs 


and increased capacity associated with digital circuits. 
[Ref. 7] 


3. Quality and Interference 


= یسک کے مت‎ e Se Se ee ____ 


The capability of digital systems to regenerate tie 
Signal and allow for multiple switching and Signal 
processing without degradation in Signal quality makes the 
digital Signal basically independent of transmission 
distance. Multiple hops from satellite to earth staticn or 
satellite to satellite are possible without accumulation of 
the ncise characteristic of analog systems. Additionally, 
digital systems are capable of operating at a signal to 
noise ratio of 20 dB to 30 dB as compared to analog systems 


requiring a much more powerful signal. (Ве. 8) 
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III. INTRODUCTION TO THE FCURIER TRANSFORM, SAMPLING 
THEO REM, 
AUTOCORRELATION FUNCTION AND THE MATCHED FILTER 
Essential to the understanding of digital signal 


modulation techniques are a few basic tools of the 
electrical engineer and the mathematician. These tools will 
be developed and elaborated on to the extent necessary to 
understand their applicability to the subject of digital 
communications. The description is not meant to be a 
detailed investigation ot the respective topics. where 
relevant, the application of the concept being described 


wili ke mentioned. 


А. FCURIER TRANS FORA 


In mathematical terms, voltages can be expressed as 
functions of time or of frequency. It is more commcn to see 
voltages represented ас а function of time as in equation 
ENIM 


V(t) = A cos (wt) | (eqn 3. 1) 


سن 
ii‏ 


amplitude 


angular frequency = 2 pi f 


frequency = 1/T 


period 


It is important tc note that both the time and £requency 
functions are representations of voltage and as such may be 
used interchangeably. The Fourier representation (v(t) ] = 


V(f) is a voitage descriptor in the frequency domain (a 


function of frequency) while v(t) is a voltage descriptor in 
the time domain (a function of time). They are different 
but equivalent and either voltage descriptor may be used, 
depending on the concept being explored, to best represent 
the voltage within context. 

The Fourier transform is of principle interest rather 
than the Fourier series since the latter is applicable only 
to periodic voltages. The Fourier transform is applicable 
to voltage pulses, random voltages and other non-periodic 
voltages. 

DEFINITION: 


ОО 
SFtv(t)] = (м) = [v(t exe jut dt (eqn 3.2) 
| -- ОО 
V(t) <-------------> У(ч) 
time domain frequency domain 


Remembering Euler's formula 
e**-jwt - cos(wt) - j sin(wt) (eqn 3.3) 


the Fcurier transform becomes 


ОО 
“7(5(5) 1 = У (м) = f vtt) [cos (wt) - j sin(vt)]  (egn 3.4) 
CO CO со ` 
У (м) = [ v(t)cos (wt) - jf v(t)ysin (wt) (eqn 3.5) 
—со — со 


Since it will be seen that digital communications deals 
priaarily with pulses of finite duration (expressed as 
period, 1118 it 15 worthwhile to examine the Fourier 
transform or a pulse cf amplitude A and duration T. 

АЕ 172 Cot < 2 
let v(t) = 


0; elsewhere 
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Figure 3.1 Square Pulse 
Figure 3.1 is a representation of v(t) or the voltage 
expressed in the time domain. The position of v(t) on the 
А 


t-axis was chosen for convenience of integration but could 


have been situated anywhere on the time line. 


CO 

Ft vty) = V(w) = [v(t)e*#-jut dt | (eqn 3.6) 
а за. T. 5 : 

V(w) = ИЕР dt + |A e**-jwt dt + (eqn 3.7) 

CO - со _ | 

Qee*s-jwt dt 2 

1/2. Tx s 

V(w) s ЈА e**-jut dt | (ez 355 
-T/L 


The integration above may te attacked either head on or 
by substituting cos(wt) - j sin(wt) for e**-jwt. The direct 
approach is illustrated due to tne relative Simplicity of 


the integrand. 


1723 
У (и) = Ја ež- jwt dt (eyn 3.9) 
зо 
T/2 
V(w) = -A/jw[e**- jwt] | (eqn 3.10) 
-T/2 
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-AÀ/jw[e**-jw1T/2 - e**Jjw1/2] (eqn 3.11) 


Ч (ч) 


Ч (м) A/jw[e** juT/2 – e**-jw«wT/2] (eqn 3. 12) 


By substituting 2pi f - w, the folloving result is obtained: 


V(f) = A/j2 pi f[e**j2 pi fT/2 - (eqn 3.13) 
e**-j2 рі ЕТ/2) 


Vif) = A/j2 pi f[e**j pi ЕТ - e**-j pi fT] 


Using Euler's formula, i-e., Sin Ө = (ежжјде-ежж—јд) /21: 


A 


V(f) = 234732 pi £[ (e**j pi fT - (eqn 3. 14( 


ежж=<ј рі ЕТ)/211 
V(f) = A/pi f{[sin( fi fT) ] 


Knowing that sin X / X = Sinc X 


V (£) АТ/рі ҒТГвіп (рі ҒТ) 1 (eqn 3. 15) 


V(f) = AT сіпс(рі ІТ) 


The sinc function iS common in digital electronics and 
plots as the product of sin(pi ЁТ) and 1/(pi ET) as in 
۲۰٢ ٦۰۴ 3:2. In Figure 3.2, V(f) in the freguency domain is 
equivalent to v(t) in the time domain. Note that as the 
pulse, 71, gets longer,  1/T gets smaller or the first zero 
crossing of the sinc function occurs at a lower and Lower 


frequency. 
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Figure 3.2 Plot of Sinc Function 


1. Amplitude and Phase Spectrum 


Recall that V(w) сап Бе represented as in equation 
ЗОНЫ The cosine term is the real part while the sine tern 
is the imaginary part. By referencing Figure 3.3, a brief 
review of the complex plane is accomplished and its 
relationship to the amplitude and phase spectrum of a given 
voltage 1s' represented. See equations 3.17 through 3.20. 


ОО CO 
V(w) = [v cos t) - JE (t) sin (wt) (eqn 3. 16) 
- ОО геа1 — COimaginary 


ІҮ (Е)| is called the amplitude spectrum or the given 
voltage. The amplitude spectrum can also be calculated by 
uSing the complex conjugate or the Fourier transform ani is 
always positive as shown in equation 3.21. The phase 
Spectrum, Ө, of the function in question is represented by 


the arctan(imaginary/real] and is illustrated in Figure 3.5. 


ED 
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Figure 3.3 Phasor Diagram 


г = [аг * b?] (eqn 3. 17) 


a + jb = |аг + b?]**.5 ежхжув (egn 3.18) 

V (f) - [real?^ * imaginary? ]**.5 e***38 (egn 3.19) 

‘1V(f) | = [rea12 + imaginary2 }**.5 (eqn 3.20) 

since |e**38| = 1 

[У (Е) | = [Y (f£) *eV*(£) ]**.5 | ыг 3. 21) 
2. Properties of the Fourier Iransform 


Several properties of the Fourier transform are 
useful in the study of digital signals. They are 
represented here without proof and without a great deal of 
detail. 

a. Linearity 
if v1(t) <-> У1(Е) апа 
if v2(t) <-> У2 (Е) then 
Flavi(t) + bv2(t) J = aV1(£) * bv2(£) 


3 1 





^i^? 





|| 


0 ےس‎ БЕ — 


Figure 3.4 Amplitude Spectrum of Square Wave 
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Figure 3.5 Phase Spectrum of Square Wave 


b.  Time-delay 
Flv (t-to) 3 = V(f) e**-jwto 
Note that the amplitude spectrum of the delayed 
versicn is the same as the amplitude spectrum oft the 
undelayed version. It is comforting to note that the 
Fourier transform of a pulse is the same tomorrow as it is 
today. 


C. Scale change 


Ftv (at) | = 1/1а| 3 (Є /а) 
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d. Frequency translation 
v(t)cos(2 pi fc t) <-> 1/2[V(ftfc) * V(f-fc)] 
v(t) is any voltage 
cos(2 pi Ес t) is a carrier wave of frequency fc 
Since understanding of this very important property 
of the Fourier transform is essential to the understanding 
of digital signal modulation it is expanded slightly here. 
As we have seen, the Fourier representation of a 
voltage pulse of amplitude A and duration T is the sinc 


function of amplitude AT as illustrated in Figure 3.6. 
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Figure 3.6 Square Wave in Time/Sinc Function in Frequency 


The translation is the product of v(t) and in this case 
cos (2 pi fc t). In the time domain this product only exists 
in the interval between -T/2 and T/2 as in Figure 3.7. The 
Significance of this translation and its relationship to the 
bandwidth of tne voltage will be discussed further in the 
section deaiing with bandwidth. 
e. Differentiation 
d v(t)/dt «-» jw V(f) 
A differentiator could be used as a clock for timing 


but would never be used in the presence of noise since the 


ЭЭ 


v(t)ces 2n f. t z(f) 


MP 
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Figure 3.7 Multiplication and Translation Diagrams 


result is exaggerated in the presence of high frequencies 
because w= 2 pi f. 
£ ° Integration 
у (+) а <-> 1/% У (5) 
An integrator could be used to reduce the effects of 
noise since w= '2 pi f is in the denominator tending to 


deeaphasize the presence of high frequency noise. 


B. THE SAMPLING THEOREM 


Essential to the understanding of digital communications 
is the sampling theorem which was first introduced by 
Nyquist in 1928 [BRef. 9] and later by Shannon in 1548 
(Ref. 10]. The sampling theorem states that any voltage can 
be uniquely represented by appropriately spaced sample 
values of the original voltage. More correctly stated, the 
Sampling theorem places limits on tne accuracy with wnicb a 
Signal can be represented. 

The implication is that an analog signal can Бе 
represented digitally or py a set of numbers, i.e., saüobles. 


A description of the sampling theorem follows. 
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A 


Given any analog voltage, v(t): as in Figure 3.8, the 
Sampling theorem says that the entire analog signal is not 
required to accurately represent the voltage but only 
Samples of it, call them vs(t). Figure 3.9 shows samples of 
a representative analog voltage. Samples can be taken of 
v(t) at every T seconds for a period of seconds. This can 
be accomplished by the use of a voltage clock, call it 
vc (t) - The sampling can be viewed grapaically as a block 
diagram representing an analog voltage multiplier as shown 
Шын Figure 3.10. To be of further use in the understanding 
of digital communications we are interested in a frequency 
description of the sample voltage, vs(t). Note that the 
System which describes the obtaining of vs(t) involves a 
voltage multiplication. It was demonstrated in the 
proceeding section that voltage multiplication amounted to 


frequency tranSlation, a property of the Fourier transforn. 
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Figure 3.8 Bepresentative Analog Voltage 


First it is necessary to find the frequency description 
of the clock, vc(t), Let vc(t) be a periodic square wave of 
height 1 and duration das in Figure 3.11. Since vc(t) is 
periodic, it can be shown 0 the Fourier series 


representing vc(t) is given by equation 3.22. 
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Figure 3.9 Samples of a Representative Voltage 


с ПС ты 


ч 


vA АУЛ ys (4) 
| | 


ус (+) 


Figure 3.10 Analog Voltage Multiplier 


Again it is noted that vs(t), the sample voltage, is the 
product of v(t), the original analog voltage, times vc(t), 
the clock voitage. In other words vs(t) = v(t) times a 
series cf cosine terms. 

оо 
vc(t) = ао +) an COS(2 pi n fc t) (eqn 3.22) 
n= 


ao and an are left unevaluated 


| 


үс (Уу 


npn 


| —————— 9 


'—(—————ÀÓ—————Á—————————————————————Á————— si m SS ара >‏ دع چ 


Figure 3. 11 A Representative Voltage Clock 


| oo 
vs(t) = у (+) Гао +.) an OS (2 pi n fc t) | (eqn 3. 23) 
n= | 


Now, assuming any form of the Fourier transform ہ٣‎ Де ДР 
as іп Figure 3.12 and since vs(t) Шоба product and by 
utilizing the frequency translation property of the Fourier 
transform, the following 1s oftained as a representation of 
the frequency spectrum of vs(t). See Figure 3.13. 

Remember that the curve nighlighted in the box in Figure 
3.13 is the Fourier represention of v(t), the original 
Signal. This original signal can now be recovered by 
filtering with an appropriate low pass filter of bandwidth 
equal to or greater than B. This low pass filter wculd only 
permit the reception of that portion of the signal which 
represents the original voltage. 

The only question left to resolve is how often to take a 
Sample. Again referring to the diagram in Figure 3.13, it 
is noted that in order to prevent any overlap of successive 


translations (aliasing) 
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Figure 3.12 Fourier Transform of v(t) 
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Figure 3.13 Fourier Transform of vs (t) 


Ғо - В >> В 

Ғо >> 2В 
or the frequency of the clock, fc (the sampling rate = 1/T) 
must Бе strictly greater than 2 times the largest 
Significant Fourier frequency component present. Since B 
will vary with different v(t), the sampling rate necessary 


to uniquely recover that v(t) also varies. 
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In summary on the sampling theorem, it can be said that 
any analog signal can be represented as certain sample 
values spaced the appropriate distance apart. The sampling 
theoren piaces linits on the accuracy of this 
representation. These sample values can then be transmitted 
from one position to another using analog to digital 
conversion and any one of a variety of modulation 
techniques. It has been shown that the Fourier frequency 
representation of the sample is not equivalent to the 
Fourier frequency representation of the source voltage; 
however, the source voltage can be recovered at the 
receiving end by filtering. A block representation of the 


system is shown in Figure 3.14. 
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Figure 3.14 Analog Voltage Multipiier and Low Pass Filter 


C. THE AUTOCORRELATICN FUNCTION 


The above concept certainly seems simple enough. If a 
signal is present and it is desired to transmit it from 
point A to point B, all that has to be done is to take 


appropriately spaced samples of the signal, encode them, 
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somehcw transmit them to point B, filter the received signal 
and the uniqueness of the original voltage has been 
reproduced. Unfortunately it isn't quite that easy. The 
reason it is not, is due to the presence of noise. The 
sources of noise will not be discussed here, however, the 
effects of noise in general terms and how a faint Signal can 
be recovered in the presence of noise will be discussed. 

The time varying descriptor of voltage and the frequency 
varying descriptor of voltage have been introduced. Both 
are precise mathematical descriptions of something (voltage) 
that is deterministic. That is, it can be described in 
mathezratical or graphical terms during any period of time. 
Such is not the case for noise corrupted voltages which are 
entirely random. Therefore other descriptors of voltages 
must be employed in the presence of noise. They are partial 
descriptors of voltage since a random signal cannot be 
described with precision. These partial descriptors are the 

1. Autocorrelation function and the 
2. Probability density function (p.d.f.) 

It is important here to note the difference between the 
Source voltage at the receiver and the Sample voltage. The 
source voltage at the receiver, vsr(t), is the digitally 
converted sample voltage, vs(t), modulated onto a carrier by 
one of the digital modulation techniques yet to be 
discussed. Understanding the autocorrelation function is 
the basis for understanding how a decision is made that 
vsr (t) is present at the receiver in the presence of noise. 
It should be remembered that the exact form of the carrier 
is kncwn at both the transmitter and receiver. 


First of all, an additive noise model is assumed where 


the signal at the receiver,  vr(t), is equal to the signal 
which is being transmitted, узг (+), (1.е., the digital 
sanples of v(t) modulated onto the carrier), plus random 


noise, vn(t) 
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vr(t) » vsr(t) * vn(t) (eqn 3.24) 
vr(t) = voltage at the receiver 

vsr (t) = signal voltage at the receiver 

vn(t) = noise voltage at the receiver 


If it is redlized that most receivers operate on the 
principal of detection of DC voltage, the problem becomes 
one of rectification of the received Signal and 
determination if the DC voltage, Characteristic of the 
transmitted signal is present. 

A common type of rectification of an analog voltage 
involves squaring the input waveforn. If the given signal 
at the receiver is vr (t) as in equation 3.24 above, squaring 
the waveform introduces the square of not only the desired 
signal but also the square of the noise term. If instead of 
squaring vr(t) and introducing or at least not eliminating 
the noise, vr(t) is multiplied by vsr(t), the results in 


equation 3.25 are obtained. 


vr (t)evsr(t) = vsr**2(t) + vsr(t) vn (t) (eqn 3.25) 


By averaging, all that remains is vsr**Z(t) since the 
average of vsr(t)vn(t) is zerc because vn(t) is random and 
the average of any random voltage is zero. Тһе DC component 
of vr(t)evsr(t) is vsr 5 Апу remaining AC component 
can be removed by a low pass filter. Graphically in block 
diagram form this receiver looks like Figure 3-15 
illustrated below. If vsr**2(t) » 0 then vsr(t) is present 
in the signal vr(t). If vr (t) is pure noise or some other 
signal is present exclusively, then the output of the 


receiver will be 0. 
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Figure 3. 15 Rectification through an AVM and LPF 
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In practice, however, the exact waveform of vsr(t) may 
not be known or it may bea time delayed or distorted 
versicn of the original Signal when it arrives at the 
receiver. If the distortion or delay is significant enough 
there will be little agreement or "correlation" in the 
receiver. | | 

The solution is to multiply the received signal by a 
series of time delayed approximations of the transmitted 
Signal. The Signal is present whenever the output of the 
Series of voltage multipliers or correlators exceeds а 
certain threshold approaching vsr**2(t). The concept is 
illustrated in Figure 3.16. 

For the concept illustrated in Figure 3.16 to work, the 


average value cf vr(t)evsr(t-p) must nave a DC component. 


vr(t)evsr(t-p) = vsr(t)vsr(t-p) + (egn 3.26) 
vn (t) vsr (t-p) 


Vr (t) evsr(t-p) 


VSI (t) VSL (t-p) + (еап 3.27) 
уп (1) узг (Е-р) 
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Figure 3.16 Correlation Device 


But the average of vn(t)vsr(t-p) = 0 since vn(t) is random and 


the results of equation 3.28 are obtained. 


vr (t) evsr(t-p) = vsr(t)vsr(t-p) (eqn 3.28) 
What 15 ОБ interest is the average value of 
vsr (t)vsr(t-p). The autocorrelation function (ACF) will be 
defined as that average value of vsr(t)vsr(t-p). The 


mathematical representation of the autocorrelation function 
for a continuous voltage is represented in equation 3.29 
while the autocorrelation function for a pulse voltage is 


presented as equation 3.30. 


: 
Rvv (p) = 1/2T | v(t) v (t-p) ät (eqn 3.29) 


RIT 
OO 


АСЕ = Суч (р) = | v(t) v(t-p) at (eqn 3.30) 
— CO 


АСЕ 


The autocorrelation function is a measure of the degree 
to which two identical Signals which are corrupted in some 
Manner are alike. The ACF of two signals which are 


identical, not distorted in any way and occurring at exactly 
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the same time, 1.е., p = 0 has a maximum value. On the 
other hand there will be very little correlation between two 
identical time variant signals when the time difference 
between them is great. In this case the autocorrelation 
function is near zero. 

A similar concept is used within the context of signal 
Comparison. This concept is the crosscorrelation which is a 
measure of the degree to which 2 different signals are 
alike. When the crosscorrelation between voltage vsr and 
random noise voltage vn is 0 there is no relationship or 
correlation. 


Again let us assume that the signal at the receiver is a 


noise disrupted version of the Signal originally 
transritted. 
vr(t) = vsr(t) + vn (t) (eqn 3.31) 


In the receiver , the time delayed version of the signal is 
applied to the incoming signal and the average is formed as 


before and Shown in Figure 3.17. 






VSr(t)vsr (4 -p)t 
—» 


уп (+) уп (+-р) 


Vr(t) >= 
vsr(t)r vn(t) 


ver (+ -p) + Vn (+-p) 
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Figure 3.17 Voltage Correlator 


44 


(vsr(t) + vn(t) ]e( vsr (t-p) * vn(t-p)] = (eqn 3. 32) 
vsr (t)vsr(t-p) * vn(t)vsr(t-p) + 
vsr(t)vn(t-p) * vn(t)vn(t-p) - 


Rvsr vsr(p) * Rvn vsr(p) * Rvsr vn(p) * Rvn vn (p) 


ii 
© 


Куп узг (р) 


Rvsr vn(p) 0; because average vn - 0 


Therefore the result is siuply Rvsr - vsr(p) * Rvn ур (р) 
or the autocorrelation function of the desired signal plus 
the autocorrelation function of the noise. Assuming the 
Shape of the autocorrelation function of both components 
(both the original signal and the noise which vary with p, 
time)is known, what is done in practice 1s to vary the value 
of p until the ratio cf Rvsr vsr(p) / Rvn vn(p) is a maximum 
or the autocorrelation function of the signal is a maximum 
while the autocorrelation of the noise iS minimum and tne 


signal is recovered. 


D. THE MATCHED FILTER 


The matched filter is a linear filter which has the 
characteristic response desired for optimum reception of the 
desired signal. In other words, we desire the response of 
the system to the linear filter to be in some way 
proportional to the autocorrelation function of the desired 
signal and in no way related to the autocorrelation function 
of the noise. This system could be illustrated as in Figure 
82213. 
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Figure 3. 18 Matched Filter Block Representation 


The question is what should the makeup or h(t),  (i.e., 
the response of the linear filter) be? The output of a 
. linear system is a convolution so the output of the linear 
filter, h(t) to an input vsr (t) would be vsr(t) * h(t), 
where * denotes convolution. 

In the previous section it was stated that the 
autocorrelation function, RVsr vsr(p), a function of p, is 
given by 

ОО 
RVSI vsr(p) = [vsr (6) уѕг (6-р) dt (ean 3.33) 
- ОО 
By Simply changing variables, equation 3.33 can be rewritten 
as a function of tine, t. 
ОО 
Rvsr vsr (t) = [уз (g)vsr (p- t) dp (eqn 3.34) 
- ОО 
Frou convoluticn theory it is known that the response of a 
linear filter as above, vsr(t) * h(t) can be expressed 
mathematically as in equation 3.35. 
GO 


усі (%) * h(t) = [vsxr(p) h(t-p) dp (eqn 3.35) 
--ОО 
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As stated earlier, the desired output of the linear 


filter is the autocorrelation function Rvsr vsr(t). Note 
the similarities between the last two equations. lets ۶۵ 
equation 3.35 a sample substitution of h(t) = vsr(-t) is 


performed, the desired results are obtained. 


ОО 


vsr(t) * h(t) = | узг(р) vsr (p-t) ар (eqn 3.36) 
- со 
- Rvsr vsr (t) 
Therefore, for every signal that is desired to be 


recovered, the matched filter will do the job very nicely if 
the response of that filter, h(t) is equal to the inverse of 


the signal that is desired to be detected. 
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IV. ANALOG TO DIGITAL CONVERSION 

Key to the understanding of digital communications is 
how the analog information 15 converted to digital 
information. As illustrated earlier, each analog voltage or 
signal can be represented by appropriately spaced sample 
values. These sample values are stiil analog in that they 
can take on any value in the range of the original analog 
signal. What is desired is to change this infinite set of 
decimal numbers to a finite set of decimal numbers. This is 
called analog to digital conversion. Common types of A to D 
conversion include Pulse Code Modulation (PCM), Differential 
Pulse Code Modulation (DPCM), Delta Modulation (DM), Pulse 
Amplitude Modulation (PAM), Pulse Duration Modulation 
(PDM),and Pulse Position Modulation (PPM). PCM, DPCM, and 
DM will ke examined in this chapter because of their 
widespread uSage and easy application to digital technology. 


A. PULSE CODE MODULATION (PCM) 


PCM 1s the most widely used A to D conversion technigue. 
It involves assigning the sample value obtained in sampling 
to one. of several quantization levels within the range of 
the voltage sampled and then representing those quantization 
levels by various binary code words. 

For exanple, the highest significant frequency component 
in human speech is 3300 Hz. In order to capture accurately 


the signal produced in speech, one would require a sampling 


rate greater than or equal to 2B or 2(3300) = 6600 
samples/sec. The telephone company uses 8000 samples/sec. 

B = 3300 Hz 

f > 28 = 6600 sanples/sec 

use £ - 8000 samples/sec to ensure no aliasing 
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If a signal is sampled at the rate of 8000 samples/sec, then 
in the transmission of those samples, assuming the samples 
are transmitted immediately and not delayed, the 


transmission rate must be 


Ts = 1/f 
Ts = 1/8000 samples/sec 
Ts = 125 microsec/sanmple 
The number of quantization levels employed to represent 
the various analog sample values 1S arbitrary. The greater 


the number of levels, the more accurate the reconstruction 
of the original signal but the more rapid the data 
transmission rate must be(in all cases there will be some 
error present after recovery). Consider again the exaaple 


of voice as illustrated in Figure 4.1. 


v(+) 





Sample every 
Т = 125 дзес 


= [= 
———————— Á— —'— 


Figure 4.1 Samples of Voltage 


Assume it is desired to represent each sample value (decimal 
number) with an 8 bit binary code word. Then the number of 
quaatization levels is given by equation 4.1. 

The spacing between quantization levels is the range of 
voltages to be represented divided by the number of 


quauütization levels. The analog to digital conversion takes 
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Figure 4.2 Analog Nature of Samples 
#juantization levels = 2**n (eqn 4.1) 


for n = 8 bits/binary code word 
#quantization levels = 2**8 


= 456 


place by determining into which quantization level the 
respective sample values fall and assigning that sample 
value the binary number associated with that quantization 
level. 

Now each sample value is assigned to an eight Lit binary 
code word. For voice transmission it has been shown that 
the Minimum transmission rate was Ts = 125 microsec/sampie. 
The bit rate is then simply equal to the number of bits per 
sample times the inverse of the transmission rate as 
illustrated in equations 4.2 and 4. 3. 

ТЕ voice was to be time multiplexed onto a single 
Channel, the bit rate for that channel would be the number 
of signals per channel times the bit rate for the most tire 


intensive signal as illustrated in equation 4.4. 


50 


bit rate - (n bits/sample) (1 sample/T sec) (eqn 4.2) 
bit rate = (8 bits/sample) (1 sample/125 microsec) 

= 64 kbps 

bit duration = 1/bit rate (eqn 4.3) 


= 1/64 kbps = 15.6 microsec/bit 


(#signals/channel) (bit rate/signal) = (eqn 4.4) 


bit rate/channel 


In other words, if 24 voice signals are to be 


multiplexed onto a single channel, the bit rate must be 


bit rate/channel 2 (eqn 4.5) 
(24 signals/channel) (64 kbps/signal) 
- 1536 Mbps 


AS can be seen, the data transmission rate increases 
Significantly, necessitating better and better hardware and 
a wider and wider bandwidth. Advantages should be readily 
apparent for A toD conversion techniques or modulation 
techniques that reduce the bit rate required. 

Recovery of the analog signal is by a process called 
Digital to Analog Conversion. All that will be said about D 
to A conversion is that it is the inverse process of À to D 
conversion and that a slight error is always introduced 
during the process due to the discrete nature of the 


quantization process during A to D and D to A conversion. 


B. DIFFERENTIAL PULSE CODE MODULATION (DPCM) 


In DPCM, what is converted to binary numbers is not the 


quantization level (decimal number) which each sample value 
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15 represented by but the successive differences between 
quantization levels. The idea is that the range of maximum 
difference values will be smaller than the range of actual 
sample values. It is therefore possible to represent that 
range of delta values with fewer bits/binary word and 
therefore fewer guantization levels and ultimately a lower 
bit rate. [Ref. 11] | 

In actual practice, in a DPCM conversion technique, гэ 
is a statistical estimate of each successive  sampie value 
which is subtracted from the actual value that is converted 
to binary code words. The result 1S the same in that the 
range of amplitudes is reduced and therefore fewer bits/word 
are required to represent the sample thereby lowering the 
data transmission rate required to transmit the signal. 
(већ. 11) 


С. DELTA MODULATION (DM) 


Delta modulation is a form of DPCM where successive 
quantization levels of the output differ by only 1 ‘bit. 
That is to say that any successive quantization level can be 
represented by varying only one bit of successive output 
binary code words. In other words a type of gray code is 
employed. This A to D technique is implemented through the 
use of a DM coder or linear delta modulator. This DM coder 
approximates a given input signal with a series cf linear 
Segments of uniform slope. A comparison is made between the 
value of this approximation andthe input signal at each 
sample increment. The sign of this difference value 1S what 
1s encoded andis used to increment the DM coder in the 
direction of the input Signal. By using the differential 
Sign value of the input signal and the incremented 
approximation from the DM coder the linear approximation 


from the linear delta modulator is said to "track" the input 
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signal. [BRef. 8] Slope overlcad of this type of modulation 
technigue occurs when the slope of the incoming signal 
exceeds the ability of the DM system to follow the source at 
the sampling rate being utilized. [Ref. 11]. 
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V. DiGITAL SIGNAL MODULATION TECHNIQUES 

Digital signal modulation techniques are the methods of 
encoding information for transmission utilizing digital 
technology. Factors affecting digital modulation include, 
but are not limited to, the physics of the method, hardware 
requirements, bandwidth considerations, power requirements, 
data transmission rates and error probabilities. It is 
these aspects which will be explored in the following 


secticns. 


A. DIGITAL MODULATION FORMATS 


Modulation is the technique by which the characteristics 
of one waveform (called the carrier) are varied or modified 
by the characteristics of another (called the source). Ihe 
carrier waveform of interest is the sinusoid. It should be 
obvious tnat the attributes of a particular sinusoid that 
differentiate it from every other  sinusoid are its 
amplitude, phase and frequency. It follows that the 
characteristics of the carrier waveform to be varied by the 
Source are its amplitude, phase, frequency . or any 
conbination of the three. This gives rise to the broad 
general formats of Amplitude Shift Keying (ASK), Phase Shift 
Keying (PSK) and Frequency Shift Keying (FSK). All other 
digital modulation techniques are variations of Or 
combinations of these basic formats. Other factors involved 
in the description of the digital modulation technique being 
empicyed inciude the number of bits being encoded at one 
tine, the employment of error correction technigues and the 


baseband (source) waveform. 
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When each bit of the baseband waveform is individually 
encoded by any of the techniques previously mentioned (ASK, 
PSK, FSK), the technique being utilized is referred to as 
binary encoding. When more than 1 bit of the source code is 
modulated onto the carrier at one time it is called block 
encoding. Block encoding allows for one of m = 2**k 
waveforms where k = number of bits. 

This paper will not go deeply into signal detection or 
demodulation techniques, however, a basic understanding of 
what is involved is necessary and again further defines the 
digital modulation format being employed. In general terms, 
signal detection is referred to as either coherent or 
noncoherent detection of the transmitted signal. Coherent 
detection, perhaps the easiest to understand, is when all 
possible waveforms of the modulated carrier waveform are 
available at the receiver and the waveform at the receiver 
15 іп phase with the transmitted carrier. Noncoherent 
detection is involved when the receiver does not have 
knowledge of the phase of the transmitted information and 
one of a number of phase estimation technigues must be 


employed for signal recovery. 


B. HARDWARE 


Although significant advances have been made in recent 
years to improve the quality of hardware associated with 
satellite communications, most components are not what would 
be considered "Off the shelf items". The hardware 
components most commonly referred to with regard to digital 
communications include the sampler, encoder, modulator, 
multiplexer and transmitter. There are variations of the 
above hardware requirements necessary for certain types of 
digital formatting, however, those exceptions will be 
addressed separately when the individual modulation 


technigues are examined. 
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Figure 5.1 Basic Hardware Component Block Diagram 


1. Sampler 


The sampler is the device or component which samples 
or extracts those single characteristic values of a 
naturally occurring analog signal which are ultimately 
encoded into a digital word by the encoder. The frequency 
at which this sampler must operate was derived in the 
discussion of the Sampling Theoren. It was determined that 
the sampler must operate at a frequency greater than 2 times 
the highest significant Fourier frequency component of the 
source voltage. Modern solid state devices capable of 
taking thousands of  samples/sec are available at modest 


costs. 
2. Encoder 


Often referred to as the A to D Converter (analog to 
digital converter), the encoder transforms the samples of 
the analog signal derived from the sampler into a digital 
format through one of the analog to digital conversion 
techniques described in the chapter on A to D conversion. 


These digital bits can be stored for later use, coded, 
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delayed or used immediately either individually or in groups 
to modify one of the characteristic qualities of the carrier 
waveform. When sanple values of the analog source signal 
аге ccnverted into digital bits of information, they are 
simply that, digital bits of information. The carrier can 
only be modified to represent the source information Бу 
interaction with another voltage. Therefore the value of 
the digital bit (binary), either 0 or 1, is used to generate 
а baseband waveform or voltage which does the actual 


modulaticn of the carrier waveform. 
a. | Common Baseband Waveforms 


It should be obvious that since it is desired to 
represent binary code words with a representative Lbaseband 
waveform what is required is two levels of voltage. There 
are two basic logic schemes for representation of the 
baseband waveforn. They are bipolar ог unipolar  iogic. 
Bipolar logic involves representing the 0's or 1's of the 


binary codeword as either +V or -V as illustrated in Figure 
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Figure 5.2 Bipolar Logic 
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The 0 or the 1 of the binary code word can be represented 
either as *V, -V or -V, +V respectively. This is a matter 
of convention but must be clearly understood in the various 
component designs in order to ensure compatibility between 
parts of the system. 

Unipolar iogic utilizes a voltage to represent 
either of the possible binary digits and the absence of 
voltage to represent the other. Two common conventions of 


unipolar logic are represented in Figures 5.3 and 5.4. 
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Figure 5.3 Unipolar Positive Logic 


Again, which form of unipolar logic might be employed in the 
A to D converter is a matter of convention. 

Variations of these two basic encoding formats 
have the advantages of ease of generation and improved 
decoding and clock recoverability. Two variations coumonly 
used in satellite communications are the Non Return to Zero 
(NEZ) and the Manchester waveforms. The NRZ waveforms is 
Simply the voltage representation which corresponds to the 
stream of bits represented in the logic scheme chosen. 
There is no transition as long as the same bit is present. 
Choosing bipolar logic as an example, a NRZ waveform can be 


described schematically as in Figure 5.5. 
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Figure 5.4 Unipolar Negative Logic 
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Figure 5.5 Non Beturn to Zero Waveform 


The Manchester waveform of the same source 
voltage is represented with a transition at the midpoint of 
the bit duration from either +V to -V or -V to +V. This 
form cf coding has the advantage of clock synchronization in 
all cases of digital encoding. Schematically, the 


Manchester code can be represented as in Figure 5.6. 
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Figure 5.6 Manchester Waveform 


Be odulator 


The modulator in a digital communications system is 
more commonly referred to as a modem 
(nodulator/demodulator). This device does the actual 
transformation of the digital information into a waveform 
that can be transmitted from one point to another. Тһе 
modulation as mentioned earlier can modify either the 
amplitude, phase or frequency of the carrier wave. In 
simple terms, the digital information derived from the 
analog source is mapped into the carrier wave. It is MES 
mapping ‘that determines the power and bandwidth 


Characteristics of the moden. 
4. Multiplexer 


Multiplexing is a common method for increasing the 
utility of a given communications channel. Variations in 
Multiplexing technigues give rise to the multiple access 


techniques employed in satellite communications. 
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a.  Prequency-Division Multiplexing (FDM) 


Frequency-Division Multiplexing is a technique 


whereby the capacity of a communications channel is 
increased by adding one signal to another. These signals 
occupy discrete nonoverlapping frequency bands. À specific 


signal is recovered by use ofa band pass fiiter for the 


frequency band desired. [Ref. 7] 
b. Time-Division Multiplexing (TDM) 


Timne-Division Multiplexing isa technique used 
to increase the capacity of a given communication channel by 
adding signals in time. That 15, specific signals or 
portions of signals are allocated specific time slots or 
portions of the carrier wave. These time allocations cannot 
overlap and the signal 15 recovered through proper 
synchronization of the recovery hardware with the 
multiplexer. [Ref. 7] 


с. Code-Division Multiplexing (CDM) 


Code-Division Multiplexing is a technique for 
increasirg the capacity of a given communications channel 
througa the assignment of a characteristic code to the 
digital signal before it is multiplexed in the time or 
frequency domain. Since the code used in multiplexing is 
Known at the demultiplexer, recovery of the digital source 
code could be accomplished through the use of a matched 


filter or correlator. [Ref. 7] 


C. BANDWIDTH 


For purposes of discusSion and for unifornity, when 
examining the various digital Signal modulation techniques 
of interest, the bandwidth will be defined as the highest 


significant Fourier frequency component of the signal in 
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question. For digital communications, it can be shown that 
by the above definition, the bandwidth of the baseband 
waveform is 1/Tb, where Tb is the bit duration. For 
example: 
1. Voice sampled at the rate of 8000 samples/sec 
2- Each sample represented by an 8 bit ccde word 
3. Rb = bit rate = (8000 samples/sec) (8 bits/sample) 
Rb = 64 kbps 
4. Tbh = bit duration = 1/Rb 
Tb = 15.6 microsec 
The Fcurier transform of the baseband waveform of duration 


15.6 microsec 1S represented in Figure 5.7. 
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Figure 5.7 Fourier Transform of Baseband Waveform 


it should бе readily apparent that as the bit rate of the 
digital baseband waveform increases, so does the bandwidth. 
It is desireable to keep the bandwidth of the signal toa 
minimum due to limitations on availability of the 
electromagnetic spectrum. In addition, interference with 
Signals which occupy adjacent frequency bands 1s reduced and 
propogation limitations of the medium and hardware 
‘Limitations such as self-inductance and capacitance are 


minimized. 
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D. SPECIFIC TECHNIQUES 


Now with some of the basics of digital signals already 
introduced it is possible to investigate some of the more 
Significant modulation formats employed in digital satellite 
communications. The organization of this section will be 
around the three general types of modulation, 1.е., phase, 
amplitude and frequency. Within each broad category, 
several techniques will be described under the subareas of 
binary encoding, where the information 1S modulated onto the 
Carrier bit by bit, and block encoding, where groups of bits 
do the modulation. Variations of these areas will be 
pointed out. The characteristics of the modulation 


technique will be introduced including, where practical, an 


analytical description, phasor representation, error 
probability, spectral efficiency, advantages and 
disadvantages. 


Phase shift keying is a generic form of signal 
modulation where the baseband waveform is used to modify or 
change the phase of the carrier. Phase, as stated earlier 
is one of tne characteristics which distinguish one sinusoid 
from another. In general, phase shift keying offers the 


advantages of power and bandwidth efficiency. 
a. Binary Phase Shift Keying (BPSK) 


AS the name suggests, Binary Phase Shift Keysing 
(BPSK) results ina mcdulated carrier waveform consisting of 
two phases or the same sinusoid. The baseband waveform is 
used to change the phase of the carrier bit by bit from one 
phase to the other. In the case of bipolar logic, the bit 
stream consists of +1. The general waveform of the carrier 


can be represented as a coSine wave of amplitude A, angular 
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frequency w and initial phase angle delta as in equation 


Sae 
ус (Е) » A cos(wt + delta) (eqn 5. 1) 
Since two phases are represented in vc(t), it is 
customary to let them differ Ly pi radians. This can be 


represented by the mcdificaticn of equation 5.1 to account 
for a pi phase shift depending on whether the bit, v(t) to 


be aodulated is either a +1 or -1 as shown in equation 5.2. 
vc (t) = À cos(wt + delta + pi/2(1-v(t))) (eqn 5.2) 


Eguation 5.2 can be expanded to resuit ina simplified form 
Showing that the modulated Carrier waveform is simply a 
product of the baseband waveform and the carrier. See 


equation 5.3 where v(t) represents the baseband waveforn. 
үс (+) = A v(t) cos(wt + delta) (eqn 5.3) 


A phasor diagram can be constructed to represent 
the modulated carrier waveform in BPSK and is shown in 
Figure 5.8. As can be seen by the phasor diagram in Figure 
5.8, the phase of the modulated waveform defends on the 
value of the baseband waveforp and differs by pi depending 
on whether the value of the modulated bit is +1. 

Remembering the section on the properties of the 
Fourier transform, it is possible to analyze the spectral 
and power efficiency of BPSK. Recall BPSK can be created by 
multiplication. Multiplication of two voltages is covered 
by the frequency translation property of the Fourier 
transform. Frequency translation results in double the 
bandwidth or spectral requirement of the translated waveform 


and half the power. In the case of a NRZ baseband waveform, 
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Figure 5.8 Phasor Diagram of BPSK 


the bandwidth of the non-translated pulse, determined from 
the largest significant Fourier frequency component is 1/Tb, 
where Тр 15 the bit duraticn. The bandwidth of the 
translated voltage (BPSK carrier) is 2/Tb. The power ina 
Sinuscid is A?/2 and none of this power capacity is lost in 
BPSK since only the phase cf the original sinusoid is 
Changed. РОГ a Manchester baseband. waveform, the 
non-translated signal has bandwidth 2/Tb and the translated 
wave has bandwidth 4/Tb. 

The advantage of BPSK is that it is relatively 
efficient in bandwidth utilization and power capacity. It 
does have the disadvantage of fairly large sidelobe 
components which contribute to interference with adjacent 
frequencies. This usually necessitates some type of 


filtering before the signal can be transmitted. 
b. Differential Binary Phase Shift Keying (DBPSK) 


The recovery of a Binary Phase Shift Keyed 
Signal requires the use of a phase coherent reference signal 
as described in tne previous sections on autocorrelation and 


Matched filters. This may not always be practical or 
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possible soa technique called Differential Binary Phase 
Shift Keying has be developed. DBPSK uses the same carrier 
type as BPSK and offers many of the same advantages and 
disadvantages. As stated, the primary difference is the 
lack of necessity for a phase coherent reference. 

The baseband waveforn, v(t), for DBPSK is 
generated by comparing the phase of the next bit to be 
transmitted to that of the previously transmitted bit. For 
example, ifa +1 is encoded onto the carrier as v(t), it 
represents a particular phase cf the carrier. If the next 
bit to be sent is aiso a *1, there is no phase change, i.e., 
the bit remains the same over two successive bit durations. 
үре however, the second bit to be sent is a - 1, this 
represents a pi phase change and is decoded as a bit change 
over successive bit durations. This can be expressed 
mathematically as in equation 5.4 where bk represents the 
bit which wili be encoded on the carrier and bk-1 represents 
the previous bit already encoded and ak represents the 


present bit. 
bk = bk-1 * ak (eqn 5.4) 


Аз сап. ре seen when Dk = -1, there is a phase 
change and the bit changes from one bit to the next 
depending on the logic type being employed. If bk ٦ 
this represents the phase remaining constant over two 
successive bit duraticns, i.e., no bit change. See Figure 
5.9 [ Ref. 8]. Recovery is effected by correlation of the 
received bit to the previously received bit. 

Differential Binary Phase Shift Keying (DBPSK) 
differs from Differentially Encoded Binary Phase Shift 
Keying (DEBPSK) in the recovery step ouly. DRPSK takes 
advantage of the modulation of phase shift in the recovery 


process while DEBPSK still utilizes a phase coherent carrier 
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Figure 5.9 Differential Binary Phase Shift Keying 
А ° 
Signal in demodulation. Both cases use the same scheme for 


differential encoding. 


с. Orthogonal Binary Phase Shift Keying (Orthogonal 
BPSK) 


In Binary Phase Shift Keying there 1S always a 
probability, with the interjection of noise on the 
transmitted signal that an individuai bit will be decoded 
incorrectiy from that which was actually transmitted. In 
some environments, this probability is so high that special 
encoding techniques must be employed in order to reduce tne 
probability that a bit or a binary code word will ke decoded 
11 error. One such technigue is Orthogonal BPSK. 
Orthogonal BPSK is a type of group encoding scheme where 
groups or blocks of bits from the original signal are 
asSigned redundant bits according to a predetermined 
assignment routine. This new expanded sequence of bits is 
then transmitted as before in BPSK. 
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For example, assume that it has been determined 
that the respective values of the original Signal will be 
represented by a k bit binary code word. It has been 
demonstrated that there exists 2**k binary code words, k 
bits long in this type of arrangement. With Orthogonal BPSK 
there also exists 2**k orthogonal sequences of bits (called 
channel symbols, bauds or chips) [Ref. 6], which represent 


the k binary code words. However, the number of bits or 
chips in the sequence is increased to n = 2**k, thereby 
providing the desired redundancy. Each of the 2**k 


sequences of n chips is constructed to be orthogoral to 
every other sequence in order to assure maximum separation. 
See Figure 5.10 [ Ref. 6]. 


=2 Orthogonal 
data word chip sequence Transmitted waveforms 
BPSK carriers with 
ohase shifts below: 
10821 Їл аны Са на тај اسیو‎ аты T 


цэ шалаа Catt) | ташилт Белі 
ар а Catt) (at) пт | 20 


оо 070 = Ca(t) 70 8 Берді 


Ел саналаа цан A ہے ےت سم ےت‎ и св 


A к ағы مس‎ Й e n m .ا‎ 


Figure 5.10 Orthogonal Binary Phase Shift Keying 


Decoding is aiso done by blocks through a bank 
ог 2*+*К correlators. Since the decoding is dore in blocks, 
the probability that the entire transmitted code word will 
be incorrectly decoded is reduced even though individual 
bits are incorrectly decoded. This is the major advantage 
of Orthogonal BPSK. Figure 5.11 [Ref. 6], Shows the 


reduction in error probability with Orthogonal BPSK 
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Figure 5.11 Probability of Error for Orthogonal BPSK 


The major disadvantage to Orthogonal BPSK is 
that. it results in a reduced data bit rate at a given 
transmission rate or the necessity for a higher transmission 
Pate af the same bit rate is to be maintained. For example, 
suppose the bit rate of a PSK signal is R bits/sec. if each 
of the k bits in the PSK Signal are represented by 2**k 


chips as is the case in Orthogonal BPSK modulation, then the 
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bit rate must be (2**k/k)R if the same data rate is to be 
maintained. Additionally, Since the necessary bit rate 
increases, the bit duration decreases. It was shown in the 
secticn on bandwidth that as bit duration decreases as the 
result of a higher bit rate, an increase in bandwidth is the 
result. Therefore, the tradeoff in decreased transmission 
error comes at the expense of bandwidth and Lit rate. 
[Ref. 6] 


d. Quadrature Phase Shift Keying (QPSK) 


Quadrature Phase Shift Keying derives its name 
from a four-phase modulation of the carrier waveform. These 
modulations are achieved by simultaneously modulating the 
inputs from two bit streams onto a singie carrier. The two 
bit sequences could be from two separate sources ог 
successive bits from a single source. To take advantage of 
the orthogonality of the sine and cosine functions, QPSK can 
be viewed as in equation 5.5 as the sum of a BPSK modulated 


sine and cosine. 


ус (+) = A vi(t) cos(wt + delta) + (eqn 5.5) 
A vz(t) sin(wt + delta) 


The signal can also be represented as in 
equation 5.6 where the phase angle, theta, of the modulated 


Carrier is shown in equation 5.7. 


vc (t) 2**.5 A cos(wt + theta + delta) (eqn 5.6) 
theta = arctan (vi(t)/v2(t)) (eqn 5.7) 
The phases possible аз a result of substitution of the 
possitie bits in birolar logic in equation 5.7 are +45 


degrees and £135 degrees. See Figure 5. 12. 
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Figure 5. 12 Modulaticn and Phases of QPSK 


Since two bits are being encoded at one time for 
transmission, the bit rate in QPSK 1S twice that of BPSK. 
The BESK bit rate is «qual to the QPSK symbol rate, (synbol 
rate being the time the guadrature bits exist on the 
Carrier). It is the symbol rate in QPSK which determines 
the bandwidth and since symbci rate ‘is equal to bit rate in 
BPSK, the bandwidth of the two modulation techniques 165 
iment ical. The great advantage of QPSK lies in the 
bandwidth utilization. Two kits of information are 
transmitted thereby effectively doubling the bit rate at no 
additional cost in bandwidth. A phasor representation of 
QPSK as illustrated in Figure 5.13. Each quadrature 
component of the modulated signal contains half the power of 
the tctal carrier or (A2/2)/2 = A2/4. 


e. Offset Quadrature Phase Shift Keying (OQPSK) 


In Quadrature Phase Shift Keying, Loth bits of 
the two baseband waveforms, which are to be modulated onto 
the carrier, change at the beginning of the respective 
Symbol time. This allows the phase of the carrier to change 


up to 180 degrees for each symbol as illustrated in the 
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Figure 5.13 Phasor Diagram of QPSK 


phasor diagram in Figure 5.13. The phase change of the 
modulated signal can be limited to 90 degrees through a 
modulation technique called Offset Quadrature Phase Shift 
Keying. This iS accomplished by delaying the application of 
the second baseband channel for Ts/2 seconds, where Ts is 
symbol duration as illustrated in Figure 5.14, #odulation 


15 then accomplished as before in QPSK. 
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Figure 5.14 Offset Quadrature Phase Shift Keying 
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The result is a QPSK modulated signal in which 
the maximum phase shift between successive bits is 90 
degrees. The bandwidth and power spectra of an OQPSK 
modulated signal are the same aS that of a QPSK modulated 
Signal. The advantage of OQESK is in the limit which is 
placed cn phase shift during encoding which simplifies the 
encoding hardware. Additionally OQPSK has spectral and 


interference advantages during decoding [Ref. 6]. 
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Figure 5. 15 Phasor Diagram of OQPSK 


f. Multiple Phase Shift Keying (MPSK) 


Another form of digital modulation is known as 
Multiple or M-ary Phase Shift Keying (MPSK). Once again, as 
150 Orthogonal Binary Phase Shift Keying, k data bits are 
transmitted. The k bits represent the word length and there 
are m = 2**k different binary words, k bits long, for this 
type of modulation. In MPSK, the binary code word is used 
to vary the phase of the carrier. There are therefore m = 
2**k different phases in MPSK. These sequences of binary 


digits need not be orthogonal as in Ortnogonal Binary Phase 
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Shift Keying. The length of the transmitted symbol is only 
k bits long as compared to  2**k bits for Orthogonal BPSK. 
The phases therefore are not all orthogonal leading to 
greater difficulty in decoding and a higher probability that 
an individual code word will be incorrectly decoded, 
especiaily as m becomes large. 

As in other types of Phase Shift Keying,  MPSK 
can be represented as in equation 5.8 where theta = phase = 


(2pi/ ği; ات ۔۔,1ء1-0‎ l. 
vc(t) = A cos (wt + theta + delta) (eqn 5.8) 


An advantage of MPSK is that as long as the 
symbol transmission rate does not increase, the carrier 
bandwidth does not increase even as the size of the binary 
code word increases. The disadvantage, as stated earlier, 
has tc de with error rates ipn decoding. Decoding requires a 
phase coherent reference of considerable Stability 
especially as m increases and respective phases become 
closer together. Practical limits for this type of 
modulation 15 m = 8 or a3 bit binary code word due to the 


complexity of the enccding and decoding hardware. 


2. Amplitude Shift Keying (ASK) 


Amplitude Shift Keying, as the name implies 1S 
modulation of the carrier through variation of its amplitude 
due tc variations in the baseband waveform. The information 
in the baseband waveform is thereby imparted to the carrier 
through this modulaticn. Since the anplitude of the carrier 
varies between successive symbols or binary code words, the 
power also varies. Recovery is by comparison of the 
transmitted power during each symbol to the possibie m - 


2**k power levels. 
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a. Multiple Amplitude Snift Keying (MASK) 


MASK is а type of group encoding where k bits 
are combined into a single waveform for transmission and 
subsequent recovery. Since each symboi or binary code word 
contains k bits, there are m - 2**k different symbols and 
consequently m different amplitudes of the carrier possible 
during the symbol duration time. An assignment scheme which 
maps the binary code words into the ш different amplitudes 
would be to simply space them Ai = A/mi volts apart, where A 
is the maximum amplitude and mi represents the decinal 
equivalent of the binary code word from i - 1 to 2**k. 

The analytic form of a MASK waveform 15 
represented in equaticn 5.9 where all the variables are the 
same a presented earlier except for Ai which is equal to one 
of i different amplitudes depending on the bit sequence to 


be transmitted. 
vc(t) = Ai cos(wt +. delta) (eqn 5.9) 


MASK 1s not very popular in satellite 


communications since it depends on a carrier of very stable 


amplitude. This is not generally practical under actual 
соп 1111 015. However, MASK does have the spectral advantage 
HS Constant bandwidth o£ 2/Ts, TS equal to symbol 


duration, even as the binary code word increases in length. 
b. QuaGcrature Amplitude Shift Keying (QASK) 


Quadrature Amplitude Shift Keying (QASK) LS a 
hybrid agigital Signal modulaticn technique. It represents 
separate amplitude modulation of the quadrature components 
of a ccmmon carrier. In that sense it is both ASK and PSK. 
Implementation ог QASK involves simultaneous application of 


M-ary Amplitude Shift Keying to the  quadratures. See 
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equation 5.10, where A1 and A2 are derived according to an 


assignment scheme as in MASK. 


ус (+) = A1 cos(wt * delta) + (eqn 5. 10) 
12 sin(wt + delta) 


By this method, the effective bit rate 15 
doubled over that of ordinary MASK. In other words, 045К 
results in the same data rate as modulating 2k bits onto the 
carrier at the same time or during the same symbol duration 
with MASK. The advantage to QASK lies in the fact that this 
increase in bit rate comes at no further increase 11 
bandwidth. The disadvantages associated with MASK are still 
present with an additional increase in the complexity of the 
decoding equipment. Decoding of QASK must be done in two 
steps. First the signal iS recovered in phase through a 
phase-coherent correlator and then each amplitude modulated 
Signal 1s recovered as before in MASK by comparison of.the 


power levels of the received signal. 


Frequency Shift Keying (FSK), is the generic term 
used to describe a number of digital signal modulation 
techniques which cause variations in the carrier frequency 
by interaction with the baseband waveforn. Decoding 1S 
generally through measurement of the frequency of the 


received signal. 


a.  Mimimum Shift Keying or rast Frequency Shift 


Keying 


Power requirement for the transmission or 
retransmission of a satellite signal 1s extremely important. 


It is desireable to limit power required to accurately 
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transmit apiece of information to the minimum possible 
consistent with allowable error rates. For this reason, 
detection of FSK signals is limited to coherent methods as 
power required increases significantly for  non-coherent 
methods. 

“SK ог ЕЕЅК are identical and represent a 
modulation technique known as continuous phase FSK. They 
get their names from the fact that "fast" indicates more 
bits per second can be transmitted in a given bandwidth than 
ordinary BPSK and "minimum" refers to the minimum modulation 
index for which orthogonal signalling occurs (Ref. 7]. 

Analytically, FFSK can be expressed as in 


equation 5.11 and equation 5.12. 


УС (1) = A cos((we + delta yw) ±) (eqn 5e T1) 


vc (t) A Cos(+deltaw t) Cos(wc t) - (eqn 5. 12) 


А sin(sxdeltaw t) sin(wc t) 


AS is noted, FFSK can be envisioned as the 
Summation cf separately modulated in-phase and quadrature 
components of tne carrier by the baseband waveforn. The 
technique is similar to that used in OQPSK. The frequency 
variation in the carrier waveform is between wc * deltaw and 
WC - deltaw. Maximum separation in phase of these two 
frequency components occurs at pi as noted in equation 5.13, 


where Tb represents the bit or symbol duration. 
(w2-w1)Tb = рі (eqn 5.13) 
This can be converted to the modulation index 


menticned above simply by converting to frequency as in 


equation 5.14. 


77 


h = (f2-£1)Tb = .5 (eqn 5.14) 


For FFSK, the frequency deviation or separation 
between frequencies of the transmitted carrier is exactly 
17/27Tb and the deviation from the carrier frequency is 1/4Tb. 
FFSK can therefore be rewritten from equation 5.12 as in 
equation 5.15 by substitution for deltaw = 2pi(1/4Tb). 


vc(t) = A cos(tpi t/2Tb) cos(2pi fc t) - (eqn 5.15) 
A sin(tpi t/2Tb) sin(2pi fc t) 


The error rate in this nodulation technique is 
identical to that of BPSK. Frequency Shift  Keying can be 
implemented to incorporate any given modulation index Бу 
sinple calculation of the frequency separation with equation 
5.14 and substitution into equation 5.15. For  M-ary 
Frequency Shift Keying, the separation between frequencies 
must be at least 1⁄/TS, where Ts is the symbol duration, in 
order to avoid carrier energy from one frequency being 
incorrectly interpreted as carrier energy from another 
frequency. This results in a modulation index of i. Figure 
5.16 represents the phasor diagram of ап FFSK modulated 


Signal. 
4. Quadrature Partial Response Signalling (QPRS) 


Quadrature Fartial Response Signalling (ОРЕ5) 
represents a specific type of the general class of digital 
Signal modulation techniques called Partial Response 
Signalling. Partial Response Signalling 1S а method in 
which a controlled amount of intersymbol interference (ISI) 
is allowed during the encoding process. Previously any 
overlap in successive Signals resulted in aliasing and an 


increased probability of error in decoding. The idea was 
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Figure 5. 16 Phasor Diagram of FFSK or MSK 


first introduced by Lender in 1963 [Ref. 12]. His concept 
was that in knowing and controlling the amount of 
interference allowed during the encoding process, 
compensations can then be made for the ISI at the receiver. 
Allowing for a limited amount of ISI makes it possible to 
transmit at rates eyual to the Nyquist rate, something that 
is not possible with many other systems [Ref. 13]. 

QPRS is implemented, as with previous guadrature 


systems, tnrough simultaneous modulation of the quadrature 


components of a common carrier. However, this time the 
modulation is accomplished with a PRS systen. The 
modulation of the respective quadrature components 


represents the impulse response of one of a number of linear 
filters to the bits cf the baseband waveform. The linear 
filter employed depends on the class of the Partiai Response 
Signalling system being used. There are 5 classes of linear 
filters commonly used in PRS and they are illustrated 
analytically in Table 1 and graphically in Figure 5.17. 
Output of a QPRS system can be represented as in 
equation 5.16, where an and bn represent the nth bit to be 


modulated and h(t - nTIs) represents the contribution of the 


да 


| ———————————————— —»——— — HQ 


li == cvm байы. meum ___-> _ EES —— 


for bite 





ОРЕ5 Linear Filter Impulse Responses 





TABLE 1 
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[ Ref. 


OO 
v< (t) z у ап h(t - n Ts) cos(wt * deita) (eqn 5. 16) 


n: -OO 
> kn h(t - n Ts) sin(wt * delta) 
n> OO 


QPRS has the advantage of rapid transmission rates 
with no increase in bandwidth and excellent error handling 
performance. The cost of such performance improvement comes 
in the form of a higher required Signal to noise (SNR) ratio 


when compared to other binary systems. 
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VI. COMPOTER SIMULATION OF DIGITAL 


Due to the nature of digital signal modulation, itis 
possikle to simulate systems through the use of the digital 
computer. This chapter 1S a description of the computer 


Simulation included in the Appendix of this thesis. 


А. | GENERAL DESCBIPTICN OF THE PROGRAM 


The program was constructed using top-down design and 
modular programming. FORTRAN was selected as the 
programming language due to its capabilities in the area of 
numerical operations and the availabiiity of mathematical 
routines to be used as modules in the progran. Testing was 
accomplished on each module as the program was developed. 
The program was designed to be used interactively in order 
to allow for instructional use. Although the program was 
nade user friendly in that error checking is accomplished on 
user inputs, care must be taken to insure instructions are 
followed correctly. 

The program consists of twenty modules, sixteen of which . 
were written by the author, three which were taken from the 
Double Precision International Mathematics and Statistics 
Library (IMSLDP) and one which was taken from a NON-IMSL 
library which resides on the IBM 3033 iocated at the Naval 
Postgraduate School. The IBM 3033 was used exclusively for 
development and testing. 

The development was accomplished using the  WATFIV 
compiler; however, the program was written to operate with 
the VS FORTRAN compiler as well. Testing and operation was 
done using the VS FORTRAN compiler. Figure 6.1 is a block 


82 


diagram of the relationship of the modules for a 


representative modulation subroutine. 


B. MAIN CONTROL MODULE 


MAIN operates as the control module for the entire 
computer simulation and accomplishes limited output. It 
introduces the program and allows the user to input various 
parameters of the digital signal modulation technique to be 
simulated and various other control functions. These inputs 
include: 

1. Modulation technique 

2. Class of QPRS system (when appropriate) 

3. Digital logic scheme 

4. Baud or symbol rate 

5. Bits/binary code word (when appropriate) 

6. Carrier freguency 

7. Number of samples to be generated 

8. Carrier max amplitude 

9. Initial phase angle 

16. Use or random number generator or no 

11. Seed for RNG (when appropriate) 

12. Number of repetitions of the simulation 

| Once the user has input the desired characteristics of 
the modulation to be simulated, MAIN calls the appropriate 
module to begin the actual calculation. The values 
necessary to calculate the time series of the signal are 
passed to the subroutines aS parameters. MAIN calls the 
reguired Subroutine the number of times the user specifies 
as repetitions. Each call to a subroutine produces the 
Discrete Fourier Transform and amplitude spectrum of the 
signal. It is these respective amplitude spectrums which 
produce the statistics in the final output after MAIN calls 


the statistics subroutine the final time. 
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Figure 6.1 Representative Block Diagram 
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C. | SUBROUTINE BPSK 


This description of SUBROUTINE  BPSK, a module that 
Simulates Binary Phase Shift Keying, will serve as a guide 
on the construction and operation of subsequent modules 
involved in the actual signal modulation. A complete 
description including its interaction with other program 
Subroutines will be included. Since other signal modulation 
subroutines interact in the same manner within the frogran, 
subseguent descriptions will concentrate on differences to 
this rasic module. 

SUBROUTINE BPSK begins with variable declarations as 
with all FORTRAN programs. Ali variables are declared to be 
either type integer or double precision. All calculations 
in the program are carried out in double precision accuracy. 

Variable initialization is next. Significant variable 
initialized at this point are TIME, the time of the first 
sample which is set to Q. STEP, the delta time between 
Samples is set to the normalized Nyquist sampling rate. 
OMEGA, the angular frequency and DELTA, the initial phase 
angle are computed in radians. NBAUD, а variable which 


keeps track of the number of symbols which have been 


modulated, is set to 1 or the first symbol is being 
modulated. ^ NBAUD remains at its .present value until total 
elapsed time exceed 1 symbol duration or  BAUDD, the 


normalized baud rate. 

Subsequently, an initial call is made to SUBROUTINE 
STREAM to receive the value of the first bit to be 
modulated. The modulation is then carried out according to 
equation 6.1. 

[he values of each sample are assigned to an array for 
storage along with the corresponding tine increment. 
Modulation continues with the first drawn bit until 1! bit or 


Symbol duration is exceeded when STREAM 1S calied to draw 
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ус (Е) = А v(t) cos(wut + delta) (eqn 6. 1) 


A = amplitude 

v(t) = bit to be modulated 
Ww = angular frequency 

t = time 


delta = initial phase angle 


another bit. This entire process is repeated until the 
number of samples specified by the user in MAIN is reacned. 

At this point, only on the first repetition of the 
Simulation, SUBROUTINE PLOT is called which gives the user 
the opportunity to plot the time series of the simulated 
Signal. Upon return from PLOT the Subroutine calls the IMSL 
routine FFTRC to generate the Discrete Fourier Transform otf 
the time series produced. This is accomplished on each 
repetition of the progran. 

Again on the first repetition only, information on the 
number. of the principie harmonic of the FFT is displayed. 
Plot is called again to give the user the opportunity to 
plot the amplitude spectrun. The values of the FFT and the 
amplitude Spectrum are computed on each Successive 
repetition of the program and the amplitude spectrum is 
added to the statistics being accumulated by a call to 
SUBROUTINE SIATS: 


D. SUBROUTINE DBPSK 


SUBROUTINE DBPSK simulated Differential Binary Phase 
Shift Keying and 1s essentially the same as SUBROUTINE BPSK. 
The difference lies in the fact that what is encoded during 
the  moduiation process is the product of the random bit 


drawn and the value of the bit previously modulated. A 
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reference bit equal to *1 is used to determine the first bit 
to be modulated. In this manner, modulation of +1 Means no 
change occurs between successive bits and the phase of the 
carrier remains the same. Modulation of a -1 indicates a 
change of bits has occurred and is indicated by a change of 
phase of the carrier between successive bit durations. The 
rest of the module remains unchanged. 


E. SUBROUTINE OBPSK 


Orthogonal Binary Phase Shift Keying is accomplished in 
this module. Again SUBROUTINE OBPSK is constructed 
essentially the same as SUBROUTINE BPSK. In Orthogonal BPSK 
the only bit streams that are modulated are n pits long 
representing each binary code word, where n = 2%*k (k = 
number of bits per word is limited to 6 bits in MAIN). Sach 
sequence of n bits is orthogonal to every other allowed 
sequence. 

The way this is accomplished in this module is first to 
draw a series of k random bits from SUBROUTINE STREAM. This 
series of k 1's and O's is then changed to its decimal 
equivalent from 0 to 2**k-1 and saved for later use. 

The program continues with a cail to SUBROUTINE ORTHO 
which generates an n Xn orthogoral matrix of +1'5 ага -1'5 
through the use of the Hadamard matrix and the Kronecker 
product [Ref. 15]. The Hadamard matrix is a 2 X 2 
Orthogonal matrix of 1's and -1's. The Kronecker product is 
the matrix which is formed when a matrix is expanded to 
twice its original size by replacing each element of the 
Hadamard matrix by the product of the Hadamard matrix and 
the original matrix. The value of the decimal equivalent to 
the k random bits is used to  identiry the row of the matrix 
and the respective columns represent the value of the binary 


digit which is modulated onto the carrier by BPSK. This 
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Sequence continues until either one bit duration 15 
exceeded, when another orthogonal bit is received from ORTHO 
or when all the orthogonal bits in a row are modulated. Then 
another sequence of k random bits is converted to its 
decimal equivalent and a new row and column of the 
orthogonal matrix is identified. Execution continues in 
this manner until all required samples have been produced. 
It shculd be noted that the bit rate for Orthogonal BPSK is 
n times the baud rate. This significantly increases the 
Signal bandwidth requirement and decreases the time 
increment between successive sampies in the simulation. The 
rest of the module 15 the same as those previously 


present ed. 


F. |SUBRCUTINE QPSK 


Quadrature Phase Shift Keying, simulated in SUBROUTINE 
ОРБК, 15 accomplished through simultaneous BPSK modulation 
of the guadrature components of the carrier. Two randon 
bits are drawn from successive calis of SUBROUTINE STREAM 
and each bit is in turn modulated onto the carrier 
components and the sur formed. This process is repeated at 
tine intervals equal to the normalized Nyquist sampling rate 
until TIME exceeds one Symbol duration, BAUDD. At this time 
two new random bits are drawn and the modulation continues 
until all required samples are produced. The rest of the 
modules 1S the same as those before. Eguation 6.2 


represents the analytical expression used to simulate 


G. SUBROUTINE OQPSK 


SUBROUTINE OQPSK is essentially the same as SUBROUTINE 


QPSK with these minor deviations. The time cr the first 
Sample is artificially initialized as TIME = .5(BAUDD) or 
1/2 the symbol duration. In this manner, the two randon 
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vc(t) = А vi(t) cos(wt + delta) + (eqn 6.2) 
A v2(t) sin (wt + delta) 


v1 (t) 
v2(t) = quadrature bit 


in-phase bit 


bits first being modulated are both known to have existed on 
the carrier for 1/2 the symbol duration. The first of thest 
two bits only remains on the carrier until .TIME = 1(BAUDD) 


when a third random bit is drawn and modulation begins with 


Ens bit. The second random bit is allowed to remain until 
TINE = 1.5(BAUDD) when a fourth random bit is drawn to 
replace it and So on. In this manner the two bits being 


modulated are offset in the time they change by 1/2 the 
svnbol duration. This necessitates keeping track of the 
number of bits modulated on both the in-phase and quadrature 
component of the carrier. Ihe remaining portion of the 
module is the same as SUBROUTINE QPSK. 


H. SUBROUTINE MPSK 


This module simulates M-ary Phase Shift Keying. This is 
accomplished by changing the phase of the carrier to any one 
Of n =2**k phases where k is the number of bits in the 
binary code word. The max length of the binary code word is 
limited to 10 in the MAIN progran. This modulation is 
accomplished by use of equation 6.3 to simulate the MPSK 
systen. 

The program first draws k random bits and converts then 
to the decimal equivalent. At this point modulation begins 
at TIME = 0 according to eguation 6.3. Modulation continues 
until one symbol duration iS exceeded when a new set of k 


random bits is drawn, conversion to decimal takes place and 
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vc(t) = A cos wt + (2рі m)/n + delta] (eqn 6. 3) 


E3 
! 


decimal equivalent of binary code word 


= 
! 


2**k; k = bits in binary code word 


modulation continues or until all desired samples have been 


produced. The program continues aS in previous subroutines. 


I. SUBROUTINE MASK 


The principles cf operation of SUBROUTINE MASK is 
Similar to SUBROUTINE MPSK which also involves a group 
encoding system. SUBROUTINE MASK allows for M-ary Amplitude 
Shift Keying. The analytic expression used in the 
Simulaticn 1S provided in equation 6.4 Ai is опе оғ 2%%қ 


equally spaced amplitudes. 
үс(%) - Ài cos(wt + delta) _ (eqn 6.4) 


Once again a set of k random bits 15 drawn . from 
SUBROUTINE STREAM. The decimal equivalent of the k bits is 
computea and used to adjust the amplitude Ai according to 


the assignment routine expressed in equation 6. 5. 


Ai = A/D (eqn 6.5) 
А = max amplitude 
поз 1 to r; n = Z**k 
This process is repeated at each symbol duration until 
all desired samples have been produced. The rest of the 


module remains unchanged. 
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4, SUBROUTINE QASK 


SUBROUTINE QASK combines the elements of the previously 
described amplitude modulation technique With the now 
familiar quadrature modulation technique. Two separate 
Streams of random numbers are generated and converted to 
their decimal forn. These numbers are used to produce two 
separate amplitudes with the same assignment scheme used in 
MASK. These amplitudes modulate the quadrature components 
of a common carrier and the output is formed by the sum of 
the quadratures. Modulation continues in the above manner 
until a symbol duraticn elapses when two new amplitudes are 
calculated. Processing terminates when all required samples 


are computed. 


Ke SUBROUTINE ASK 

Minimum Shift Keying is frequency shift keying in which 
the modulation index is 1/2. The modulation index is 
represented in equaticn 6.6 


h = „5 (deltaf) Tb (eqn 6.6) 


deltaf = frequency separation 
Tb = bit duration 


The frequency deviation frcm the carrier is .5(deltaf) - 
174. Converting this frequency deviation to radians 
yields pi/2Tb. This angular frequency deviation from the 
central carrier is either +pi/1Tb depending on the value of 
the bit to be modulated. 

The simulation in SUBROUTINE MSK 1s accomplished ру 
drawing a random bit and generating the signal according to 


equation 6.7. 
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vC (t) = À cos( (Ww +pi/2Tb)t + delta] (ечп 6.7) 


Once again modulation continues until one bit duration 
is elapsed when another bit is drawn. Modulation stops when 
all required samples are produced. The rest of the 


Subroutine remains unchanged. 


L. SUBROUTINE МЕ5К 


.This subroutine modulates a signal simulating  M-ary 


Frequency Shift Keying. The separation between adjacent 
freyuencies is estabiished as 1/Ts, where Ts is the symbol 
duration. This amounts to a moduiation index egual to 1. 


Initially the ertire range of frequency deviation from the 
Carrier is calculated as (n-1)Ts. The mean frequency 
deviation is then calcuiated. In other words the range of 
frequencies of the modulated signal is the frequency of the 
central carrier +1/2 the magnitude of the entire range of 
frequency deviation. The minimum frequency is used as the 
base fcr computing the frequency to be used to modulate the 
signal during each respective symboi duration. | 

At this point in the subroutine a series of  k random 
bits (k<10) 1S drawn end converted to decimal.  Tbis decimal 
number is multiplied by the frequency Separation and added 
to the minimum frequency and converted to radiars. The 
analytical expressicn of the signal produced with this 


subroutine is illustrated in eguation 6.8. 


vc(t) = À cos[ (q + 2pi mf)t + delta] (eqn 6.8) 
mf - moduiation frequency 
This modulation continues until a symbol duration 


elapses when a new frequency deviation or mf is calculated 
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or mcdulation stops do to completion of all required 


samples. 


N. SUBROUTINE ОРВ5 


This mcdule simulates 5 different classes of Quadrature 
Partial Response Signalling systems. The technique employed 
involves the summation, at each time increment, of all 
responses to the linear fiiter representing the class. This 
is accomplished for all of the bits being modulated during 
the duration of the simulation. 

Initally, two arrays of random bits are generated 
representing those bits which would be modulated onto the 
quadrature components of the carrier during the length of 
time the signal is to be simulated. The the impulse 
response to each of these n bits is calculated for the time 
of the respective sample. These impulse responses represent 
the modulation for the bit in question onto the respective 
porticns of the carrier and the sum is formed according to 


equation 6.9. 
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Ус (+) = УА h(t - п15) cos (wt + delta) + (eqn 6.9) 
CO n:-Oo 

У А h(t - nTs) sin (wt * delta) 

--ОО 

h(t - nfs) = impulse response of nth bit at time t 


The time of the first sample is initialized at 6(BAUDD) 
in order to ensure that when the first sample is taken, 
contributions from bits modulated at time 0 are also summed 
do to the overlapping nature of the QPR signals. The value 
of the last sample is also formed by the sum of the 
quadrature components existing at 6(BAUDD) after it 15 
produced. The rest of the program operation remains the 


same as previous modules. 


ELS 


|. SUBROUTINE STREAS 


This subroutine interacts with the IMSL random number 
generator SUBROUTINE GGUBFS to produce a random nunber 
between 0 and 1 and make assignment on the basis of the 
value of this random number to random bits according to the 
digital logic scheme specified in the MAIN program. [It also 
enables the user to manually insert binary digits if that 


option was specified previously in MAIN. 


O. SUBROUTINE ORTHO 


СБТНО produces an n Xn orthogonal matrix from the 
Hadamard matrix by forming successive Kronecker products 
fret. m5]. It also selects the appropriate row and coluan 
of the orthogonal pit to be soduiated and passes this bit 
back to SUBROUTINE OFPSX. Foint of entry to SUBROUTINE 
ORTHO 1S controlled by a flag set in OBPSK. 


P. SUBROUTINE PLOT 


PIOI 1S an interactive module that allows the user to 
determine whether or not a graph of the output of program 
calculations is to be produced. The user is also able to 
selectively vary certain parameters associated with the 
plot. PLOT calls SUBROUTINE UTPLOT which actually produces 
the graph and performs the output. JTPLOT 1s a NON IT 


library routine from the Naval Postgraduate School. 


Q. SUBROUTINE STATS 


The statistical calculatiors associated with the 
amplitude spectrums of the various functions generated on 
successive repetitions of the program are computed іп 
SUBROUTINE STATS. Upon coupletion of each repetition, each 


module calls STATS where the Sus, sum of the squares, sum of 


the cubes and sum of the quartics of each element of the 
авріі обе spectrum of the Signal produced are computed and 
stored. when the last repetition of the prograz is 
complete, БАГИ instructs STATS to compute and output the 
Гїлай statistics, і.е., mean, Variance, skewness and 
kurtosis оі each component of tke amplitude spectrin 
according to the estimations contained in [Ref. 16]. 
Additionally the wean variance and variance of the variance 
is calculated and  cutput. Point ог entrance to this 
subroutine 15 controlled by flags set inthe individual 


eodulaticn subroutines or in MAIN. 
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VII. STATISTICS OF IHE PAST POURIER 


A. DESCRIPTION OF THE PROBLEM 


As stated in Chapter II, if it can be shown that the 
statistics of the FFT can be somehow linked to a particular 
digital signal modulation technigue, then the modulation 
technique can be identified on the basis of those statistics 
alone. Although it 15 beyond the scope of this thesis to 
actually derive those relationships if they exist, it seen 
prudent to attempt to determine if there is statistical 
differences between the components of the FFT's as they are 
derived in the enclosed computer simulations. 

The statistic chosen to do the hypothesis testing is the 
F-test since the true mean and variance of the distribution 
need not be known [Ref. 17]. In addition, the information 
necessary to calculate the F-statistic 15 readily 
accumulated in  SUBROUTINE STATS. Calculation of the 
F-statistic necessitated the writing of a small computer 
program to be used for that purpose once the output from the 


main frogram was generated and reformatted. 


B. ANALYSIS—OF-VARIANCE 


If there were no differences attributable to the 
modulation technique, then it would be reasonable to assume 
that for a certain set of characteristics, the components of 
successive FFT's would be the same. In other words, if a 
Signal was nodulated by BPSK and a statistically significant 
number of FFT's were generated of the signal, then the 
Statistics of those rFT's would be assumed to be related. 
On the other hand, if it can be shown that the statistics 


are not related (i.¢€., not from the same distribution), then 
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significance can be placed in the variation among modulation 


techniques. 


1. The F-Distribution 


The assumptions necessary to use the F-distribution 

are: 

1. Normal distribution 

2. Random samples 

3. Independent samples 
These assumptions are not too difficult to intuitively 
accept in tne model that will be proposed. Once again it 
would be expected that the FFI's of successive identically 
modulated signals to be related. For the simulations 
conducted as part of this test, they were all generated 
using bits from a random number generator for which it can 
be shcwn that each sequence of bits passes statistical tests 
for  randomness and independence. The assumption of 
normality could also be argued on the basis of the central 
limit theoren. 

The vaiue of the F-statistic with which comparison 
will be made is 1.51 for 14 degrees of freedom іп the 
humerator, an infinite number of degrees of freedom in the 
denominator and a 90% confidence region. How these values 
were oktained will be detailed under the design of the 


experiment. 


2. Design of the Experiment 
The experimental data used in the computation of the 
Statistics was chosen to minimize the random contributions 
of variables other than the modulation technicue. The 
computer simulations included in this thesis were used to 
generate the statistics and a separate program also included 
in the appendix was used to calculate the F-statistic. The 


variables in the experiment include: 
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1. Modulation technique 

2. Logic type 

3. Baud rate 

4. Bits per binary code word 

5. Carrier frequency 

6. Carrier amplitude 

7. Initial phase angle 

8. Number of samples generated 

9. Time between samples 

10. Seed for random number generator 

11. Number of repetitions or trials 

Fifteen different modulation techniques were 
compared. in all cases bipolar logic was simulated, the 
baud rate was held constant at 1200 baud and the maximum 
carrier amplitude was established at 1 volt. Additionaily 
the phase angle of ail simulations was 0 degrees and the 
seed for the random number generator was 1 thereby ensuring 
the same random sequence of bits. The number of bits per 
binary code word did vary among the modulation techniques. 
When the modulation technique was M-ary, the bits per code 
word was always 3 so it is possible to infer that bits per 
code word 1s a function of the modulation technique. The 
time between samples was the normalized Nyguist sampling 
rate. This did vary between modulation techniques but was 
necessary to derive an accurate FFT. Carrier frequency also 
varied between Simulations but was always twice the lowest 
Carrier frequency recommended in the computer Simulation. 
In most cases this was 2400 Hz, or twice the baud rate. 
Finally each sinulation was repeated 100 times and the 
Statistics of the FFT's based on those 100 repetitions. The 
number of samples produced was always 64 so there are 33 
components of the FFT. 
What this amounts to can be illustrated through an 


Analysis-of-Variance table as shown in Table 2 [Ref. 17]. 
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TABLE 2 
ANALYSIS-OF-VARIANCE TABLE 


n Observations in Each of r Groups 


[Observations | зөв | Mean 
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In the case of the simulation described above, this amounts 
to 100 observations (repetitions or triais) from each of 15 
groups (modulation techniques) for each of the 33 components 
of the FFT. The degrees of freedom in the  numerator is 
equal to (15-1) or 14 while the degrees of freedcm in the 
denominator is equal to 15(100-1) or 1485. Table values fcr 
the F-distribution use infinity as the degree of freedom for 
values greater than 120. 

An F-test was also performed to determine if there 
15 a relationship among the mean variance, mean skewness or 
mean kurtosis of the 15 different modulation techniques. 
Again the assumption is made that these statistics wouid be 
normaliy distributed from modulation technigues which had 
the same or statistically similar components of their FFT's. 
In this case there were again 15 nodulation techniques to 
compare (14 degrees of freedon) but only 33 observations 
(the variance, skewness or kurtosis of the respective 
elements of the FFT). This number ot observations yields 


15(33-1) = 480 degrees of freedom in the denominator. 
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3. Hypothesis and Hypothesis Testing 





The hypothesis posed for the model is that the means 
of the individual ccmponents of the  FFT's from the 15 
modulation techniques are fron the same distribution. 
Therefore, the respective means must Бе equal. Also tested 
is that the mean variance, mean skewness and mean kurtosis 
of all components of the  FFT's of each modulation technique 
are equal. These hypothesis are illustrated in equations 
7.1 through equation 7.4. Each hypothesis is tested and 
compared to the F-distribution selecting a rejection region 
or level of significance of .9. This equates to an 
F-statistic of 1.51 for 14 degrees of freedom in the 
numerator and an infinite number of degrees of freedom in 


the denominator. 
4. Results 


The means of the 33 ccmponents of the FFT for each 
of the 15 modulation technigueS were compared using the 
F-distribution. The results are shown in Tabie 3. Table 4 
Shows the F-statistics associated with the comparison of the 
mean of the variance, mean of the skewness and mean of the 
kurtosis for the respective groups. In addition a complete 
summation of the results of the F-test are included in the 


appendix. 


C. CCNCLUSION 


The results of the F-test indicate that the null 
hypothesis uust be rejected at the .9 Significance level and 
the aiternate hypothesis accepted that the means are not 
equal for FFT components 9-14. In addition the F-statistic 
for the comparison of mean variance, mean skewness and mean 


kurtosis all fall in the rejection area. 
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Ho: XBAR11 = XBAR21 = ө ө • = ХВА81] (eqn 7.1) 


H1: XBAR11 «4. ХВАВ21 #4 e e ө Ұ ХВАВ13 


P. 
! 


ith modulation technique 
jth component of the FFT 


Ho: MVAR1 = MVAR2 MVAR15 (eqn 7.2) 


H1: MVAR1 * HVAR2 * ee # МУДЕ 15 


HO: MSKEW1 = MSKEW2 = e è • = MSKEM15 (eqn 7.3) 
Н1: М5КЕН1 # МЅКЕН2 # e ө ө # MSKEW15 
Ho: MKUR1 = MKOR2 = ө ө e = MKURI5 (eqn 7.4) 


H1: МКОКТ + МКОК2 4 e e # MKURÍ5 


This indicates that there is a difference between the 
Statistics of of the FFT's of the respective modulation 
technigue. Since all the parameters used in the generation 
of these statistics were held essentially constant among the 
tiaals, it can be assumed that these differences are due to 
the modulation technique employed. It has not been 
established yet what those differences may be. This is an 
area where future research and study is warranted. 

The results of this test would seem to point to those 
components of the FFT which offer the best opportunity to 
develcp those relationships or differences. FFT components 
9-14 have oktvious differences; however, FFT components 20, 
24 and 25 also have large F-statistics but do not fall in 


the rejection area. These components may also be 


10 1 


OEELA و سے ہے وھ ہے مک ہے‎ O Oo ety MIR etn ps GN pee, EN, SE GER —— GEE pe, АШЫҒЫ a. — — SR бару ОиВ те: 


TABLE 3 
F-STATISTICS OF FFT COMPONENTS 


Lc m нанели کس‎ RES ON Qm S am, dame RR UNS سے‎ шинж шиль amas 


| 

| 

| 

| 

| FFT AMONG WITHIN TOTAL 

| | CCMPCNENT GROUP SUM GROUP SUM OF SQUARES  F-STAT 
P" 1 02121 17.69 17.81 0. 729 
| 2 0.112 18.48 18.59 0.642 
| 3 0.108 20.65 20.76 0.557 
| Т 0.135 25.40 25.54 0.562 
| 5 0.172 32.74 32.91 0.558 
| 6 0.254 41.72 41.98 0.647 
| 7 0.422 50.10 50.52 0.893 
| 8 0.962 82.83 83.79 1.231 
| 9 7.591 390.2 397.8 2.063 
| 10 49.09 1681.0 1730.1 3.097 
11 53.01 1732.1 1785.1 3.246 
| 12 61. 56 2071.7 2133.3 3.152 
| 13 45.73 1832.8 1878.6 2.646 
| 14 3.864 250.5 254.4 1.636 
| 15 0.539 98.13 98.67 0.582 
| 16 0.366 110.1 5۔110‎ 0. 353 
| 17 0.504 151.0 151.5 0.354 
| 18 0.943 140.3 1:81 ۔‎ 2 0.7.13 
| 19 1.299 158.5 159.8 0.370 
| 20 2.671 200.7 203.4 1.912 
| 21 1.853 226.3 228.2 0.868 
| 22 22 232.1 234.3 1.041 
| 23 1.712 197. 6 199.3 0.919 
| 24 2.512 18 1. 9 184.4 1.476 
| 25 2.144 17 1.3 173. 5 1. 327 
| 26 1.257 143.7 145.0 0.928 
| 27 1.233 113.9 115.1 1.149 
| 

| 











ма” NEL T "71 
| TABLE 3 | 
| F-STATISTICS OF FFT COMPONENTS (con't) | 
| 28 0.858 115.8 116.7 0.785 | 
| 29 0.547 63.72 64.27 0.911 | 
| 30 0.356 5527? 56.13 0.677 | 
| 31 0.183 41.41 51.59 0.469 | 
| 32 0mO 51 13.79 13.84 0.389 | 
| 33 0.031 10.92 10.96 0.34 | 
E Е. d 
| л. ээ | 
| TABLE 4 | 
| F-STATISTICS OF THE MEAN VARIANCES, | 

MEAN SKEWNESS AND MEAN KIRTOSIS | 
| | 
| AMONG Wi THIN TOTAL Í 
| GROUP SUM GROUP SUM OF SQUARES F-STAT | 
| VARIANCE 3.144 25 3.755 E6 4.070 Еб 2.870 | 
| SKEWNESS | 5.307 E10 7.031 E11 fee oe me ie # 22588 | 
| KURTOSIS 8.280 E11 1.262 E13 0 0 | 
—— eee 17 | 


In addition, it should be noted 


that the statistics of the FFT associated with each 


interesting to examine. 
respective modulation technique also display some striking 


differences. These statistics may prove in some way to be a 


fingerprint of the modulation techniques themselves. 
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APPENDIX E 

RESULTS OF F-TEST 
F-STATESTTCO 77777 1 COMPONENT 
SUM OF SUMS = 0.2091487882294994D 02 
SUM OF SUMS**2 = 0.1810655469453404D 02 
SUM OF SQUARE OF SUMS = 0.4128097510679048D 02 
AMCNG GROUP SUM OF SQUARES = 0.1211883136154515D 00 
WITHIN GROUP SUM OF SQUARES - 0.1769374494346614р 02 
TOTAL SUM OF SQUARES = 0.1781493325708159D 02 
DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.148500000C000000D 04 
F-STATISTIC = 0.7265063214387600D 00 


F-STATISTICS* FOR THE 2 COMPONENT 

SUM OF SUMS = 0.2305792717805020D 02 

SUM OF SUMS**2 = 0.1894701234015291р 02 

SUM OF SQUARE OF SUMS = 0.4663825126655468D 02 

AMCNG GROUP SUM OF SQUARES = 0.1119371755000360D 00 

WITHIN GROUP SUM OF SQUARES = 0. 1848060298274 8736D 02 

TOTAL SUM OF SQUARES = 0.1859256709298739D 02 

DEGREES OF FREEDOM IN NUMERATOR - 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 
Р-ЗТАЧТЬТТС = 0.64247464650 147370 00 


F-STATISTICS FOR THE 3 COMPONENT 

SUM OF SUMS = 0.2693734535827635D 02 

SUM OF SUMS**2 = 0.2124674921974186D 02 

SUM OF SQUARE OF SUMS = 0.5921566381358867D 02 

AMCNG GROUP SUM OF SQUARES = 0. 108409588 1685 183р 00 

WITHIN GROUP SUM OF SQUARES = 0.2065459258160598D 02 

TOTAL SUM OF SQUARES = 0.2076300216977449D 02 

DEGREES OF FREEDOM IN НОМЕВАТОЕ < 0.14000000000000000 02 
DEGREES OF FREEDOM IN DENOMINATOR < 0.14850000000000090 04 
F-STATISTIC = 0.5567362242775778D 00 
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Е-ЭГМЕР5ТІС5 GOR THE 4 COMPONENT 

SUM OF SUMS = 0.3175688631627065D 02 

SUM OF SUMS**2 - 0.2621104542300214D 02 

SUM OF SQUARE OF SUMS = 0.8070513093726266D 02 

AMONG GROUP SUM OF SQUARES = 0.1347180903696015D 00 

WITHIN GROUP SUM OF SQUARES = 0.2540399411362951D 02 

TOTAL SUM OF SQUARES = 0. 2553871220399911D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0. 14850000000000 00D 04 


F-STATISTIC = 0.56249974850418690 00 

КӘДЕЛІ О ТРС5 FOR THE 5 COMPONENT 

SUM OF 5045 = 0.3669104349531116D 02 

SUM CF SUMS**2 - 0.3380633740317695D 02 

SUM OF SQUARE OF SUMS - 0.1069682473356158р 03 

AMCNG GROUP SUM OF SQUARES = 0.1721940248396148D 0 
AITHIN GROUP SUM OF SQUARES = 0.32/3665492982079D 02 
TOTAL SUM OF SQUARES = 0.3290884895466041D 02 

DEGREES Or FREEDOM IN NUMERATOR = 0.1400009000000000D 02 
DEGREES CF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.5579331866788867D 00 


H STATISTICS FOR THE 6 COMPONENT 

SUM OF SUMS = 0.41600978 789484310 02 

SUM OF SUMS**zZ = 0.43132366420200940 02 

SUM CF SQUARE OF SUMS - 0.1408215358007360D 03 

AMCNG GROUP SUM OF SQUARES = 0.2584544005119440D 00 

WITHIN GROUP SUM OF SQUARES = 0.41724151062193580 02 

ТСТАІ SUM OF SQUARES = 0.4197860546270552D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 92 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.6468757561623442D 00 
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F-STATISTICS FOR THE 7 COMPONENT 

SUM CF SUMS = 0.4810387092348111D 02 

SUM OF SUMS**2 = 0.52064093 81611744D 02 

SUM OF SQUARE OF SUMS = 0.1964653659903027D 03 

AMONG GROUP SUM OF SQUARES = 0.4219987280210724D 00 

WITHIN GROUP SUM OF SQUARES = 0.5009944015621441D 02 

TOTAL SUM OF SQUARES = 0.5052143888423548D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.8934632362547178р 00 


Е-ЕТАТІ5ТІС5 FOR THE 8 COMPONENT 

SUM OF SUMS = 0.5623255941290002D 02 

SUM CF SUMS**2 = 0.8590285622098 203D 02 

SUM CF SQUARE OF SUMS = 0.3070018365343552D 03 

AMCNG GROUP SUM OF SQUARES = 0.9619512065933320D 00 

WITHIN GROUP SUM OF SQUARES = 0.8283283785563848D 02 

TOTAL SUM OF SQUARES = 0.8379478906223181D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.1231824736914027D 01 


F-STATISTICS FOR THE 9 COMPONENT 
SUM OF SUMS = 0.9931605295569187D 02 

SUM OF SUMS**2 = 0.4043531056175881D 03 

SUE OF SQUARE OF SUMS = 0.1416642417283150D 04 

AMCNG GROUP SUM OF SQUARES = 0.7590638589699638D 01 

WITHIN GROUP SUM OF SQUARES = 0.3901866814447566D 03 

TOTAL SUM OF SQUARES = 0.3977773200344562D 03 

DEGREES OF FREEDOM IN NUNERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0. 148500000G000000D 04 
F-STATISTIC = 0.2063499133280512D 01 
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F 500000 500 کت‎ FOR THE 10 COMPONENT 

SUM OF SUMS = 0.2029942628781924D 03 

SUM OF SUMS**2 = 0.1757564757819924D 04 

SUM OF SQUARE OF SUMS = 0.7655692250451960D 04 

AMCNG GROUP SUM OF SQUARES = 0.4908580866354581D 02 

WITHIN GROUP SUM OF SQUARES = 0.1681007835315404D Ow 

TOTAL SUM OF SQUARES = 0.1730093643978950D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.14850000000000000 04 
F-STATISTIC = 0.30$7309684192641D 01 


D- STATISTICS FOR THE 11 COMPONENT 

SUM OF SUMS = 0.2004773302247246D 03 

SUM CF SUMS**2 - 0.1811875077561617D 04 

SUM OF SQUARE OF SUMS = 0.7980439863108266D 04 

AMCNG GROUP SUM OF SQUARES = 0.5301029200839379D 02 

WITHIN GROUP SUM OF SQUARES = 0.1732070678930534D 04 

TOTAL SUM OF SQUARES = 0.1785080970938928D 04 

DEGREES OF FREEDOM IN NUMERATOR - 0.1400000000000600D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.32463325375329120 01 


A STATISTICS POR THE 12 COMPONENT 

SUM OF SUMS = 0.2210834652773109D 03 

SUM CF SUMS**2 = 0.2165859765218944D 04 

SUM OF SQUARE OF SUMS = 0.9414279932820438D 04 

AMCNG GROUP SUM OF SQUARES = 0.6155753358218842D 02 

WIIHIN GROUP SUM OF SQUARES - 0.2071716965890740D 04 

TOTAL SUM OF SQUARES = 0.2133274499472928D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0.14000000000090000 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.14850000000900000 04 
F-STATISTIC = 0.3151731454585565D 01 
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F-STATISTACS FOR THE 13 COBPONENT 

SUM OF SUMS = 0.1979583434471791D 03 

SUM OF SUMS**2 + 0.19047010 753471600 04 

SUM OF SQUARE OF SUMS < 0.7 1854466 728408620 04 

AMCNG GROUP SUM OF SQUARES = 0.4572946290150773D 02 

WITHIN GROUP SUM OF SQUARES = 0. 18328466086 187510 04 

TOTAL SUM OF SQUARES = 0.1878576071520259D 04 

DEGREES OF FREEDOM IN NUMERATOR = 0. 1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D OY 
F-STATISTIC = 0.26464786 714598640 01 


F-STATISTICS FOR THE 14 COMPONENT 

SUM CF SUMS = 0.92277362 726 146190 02 

SUM OF SUMS**2 + 0.2600827156 8862540 03 

SUM OF SQUARE OF SUMS - 0.9540726819081072D 03 

AMCNG GROUP SUM OF SQUARES = 0.38063985704619236D 01 

WITHIN GROUP SUM OF SQUARES = 0.250541992869 5444) 03 

TOTAL SUM OF SQUARES = 0.2544059785741636D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000006000000D 04 
- F-STATISTIC = 0.1635887377498235D 01 


F-STATISTICS FOR THE 15 COMPONENT 

SUM OF SUMS = 0.6880926336210145D 02 

SUM CF SUMS**2 - 0.1018300362534556D 03 

SUM OF SQUARE OF SUMS = 0.3695062863415098D 03 

AMCNG GROUP SUM OF SQUARES = 0.5385863 804584 070D 00 

WITHIN GhOUP SUM OF SQUARES = 0.9813497339004056D 02 

TOTAL SUM OF SQUARES = 0.9867355977049896D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000р 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.14850000000002900D 04 
F-STATISTIC = 0.5821433971075596D 00 
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F-STATISTICS FOR THE 16 COMPONENT 
SUM OF SUMS = 0.6961387326450431р 02 

SUM OF SUMS**2 = 0.1137157033969463D 03 

SUM OF SQUARE OF SUMS = 0.3597053597878584D 03 

AMCNG GROUP SUM OF SQUARES = 0.3663260306209388D 00 

WITHIN GROUP SUM OF SQUARES = 0.1101186497990677D 03 

TOTAL SUM OF SQUARES = 0.1104849758296886D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC - 0.3528623485827816D 00 


F-STATISTICS FOR THE 17 COMPONENT 

SUM GF SUMS = 0.7052854165691613D 02 

SUM OF SUMS¥*2 = 0.1548042004914454D 03 

SUM OF SQUARE OF SUMS - 0:3820324148469269D 03 

AMCNG GROUP SUM OF SQUARES = 0.5041406896350342D 00 
WITHIN GROUP SUM OF SQUARES = 0.15098387634297610 03 
TOTAL SUM OF SQUARES = 0.1514880170326111D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.3541763825767676D 00 | 


Босо FOR THE 18 COMPONENT 

SUM OF SUMS = 0.7950452893326030D 02 

SUM CF SUMS**2 + 0.14541276561222910 03 

SUM OF SQUARE OF SUMS = 0.5 157161626920556р 03 

AMCNG GROUP SUM OF SQUARES = 0.9431815463208067D 00 

WITHIN GROUP SUM OF SQUARES = 0.1402556039853086D 03 

TOTAL SUM OF SQUARES = 0.14 11987855316294р 03 

DEGREES OF FREEDOM IN NUMEKATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.7133020797581566D 00 
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F-STATISTICS FOR THE 19 7٤٦ 

SUM OF SUMS - 0.8206200282896035D 02 

SUM OF SUMS**2 - 0.1642503769158388D 03 

SUM OF SQUARE OF SUMS = 0.5788746875066006D 03 

AMCNG GROUP SUM OF SQUARES = 0.1299298669532475D 1 

WITHIN GROUP SUM OF SQUARES = 0. 1584616300407728D 03 

TOTAL SUM OF SQUARES = 0.1597609287103052D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 


F-STATISTIC = 0.8697276809711277D 00 


F-STATISTICS FOR THE 20 COMPONENT 

SUM OF SUMS = 0.9520048302882186D 02 

SUM OF SUMS**2 = 0.2094086558414699D 03 

SUM OF SQUARE OF-SUMS = 0.8713219218815827D 03 

AMCNG GROUP SUM OF SQUARES = 0. 2671131239535162D 01 
WITHIN GROUP SUM OF SQUARES = 0.2006954366226541D 03 
TOTAL SUM OF SQUARES = 0.2033665678621893D 03 ` 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000099D 02 
DEGREES OF FREEDOM IN DENOMINATOR з 0.14850000000000000 04 
F-STATISTIC = 0.14117446377815830 01 | 
F-STATISTICS FOR THE 21 СОМРОМЕНТ 


SUM OF SUMS - 0.8994139952150146D 02 

SUM OF SUMS**2 - 0.2335862791765991D 03 

SUM OF SQUARE OF SUMS = 0.7245755635510632D 03 

AMCNG GROUP SUM OF SQUARES = 0.1852785403586404D 01 

WITHIN GROUP SUM OF SQUARES = 0.2263405235410884D 03 

TOTAL SUM OF SQUARES = 0.2281933089446748D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR - 0.148500000C000000LD 04 
F-STATISTIC = 0. 8682828488687504D 00 
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P- STATISTICS FOR THE 22 COMPONENT 

SUM OF SUMS = 0.1035648300735959D 03 

SUM OF SUMS**2 - 0.2014810868436007D 03 

SUM OF SQUARE OF SUMS = 0.9427341493765320D 03 

AHCNG GROUP SUM OF SQUARES = 0.2276892141650119D 01 

WITHIN GROUP SUM OF SQUARES = 0.2320537453498354D 03 

TOTAL SUM OF SQUARES = 0.2343306374914855D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = Q0.1485000000000000D 04 
F-STATISTIC = 0.1040764077320930D 01 


F—STATISTICS FOR THE 23 COMPONENT 

SUM OF SUMS + 0.9031234717154 7460 02 

SUM OF SUMS**2 = 0.2047768116906234р 03 

SUM OF SQUARE OF SUMS = 0.7149379755971323D 03 

AMCNG GROUP SUM OF SQUARES = 0.17118330548819120 01 

WITHIN GROUP SUM OF SQUARES = 0.1976274319346€521D 03 

TOTAL SUM OF SQUARES = 0.1993392649895340D 03 

DEGREES OF FREEDOM IN NUMERATOR - 0.1400000000000000D 02. 
DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000р 04 
F-STATISTIC = 0.91878 22552243536D 00 


Е-5ТАТІ5ТІС5 FOR THE 24 COMPONENT 

SUM OF SUMS - 0.9171706807335522D 02 

SUM OF SUMS**2 - 0.1900438745495707D 03 

SUM CF SQUARE OF SUMS = 0.8139562480340389р 03 

AMCNG GROUP SUM OF SQUARES = 0. 2531548763025405D 01 

WITHZN GROUP SUM OF SQUARES = 0.1819043120692303D 03 

TOTAL SUM ОЕ SQUARES = 0.1844358008322557D 03 

DEGREES OF FREEDOM IN NUMERATOR - 0.1400000000000000D 92 
DEGREES CF FREEDOM IN DENOMINATOR = 0.1748500000C00000UD 0% 
F-STATISTIC - 0.1476188171339999Dp 01 
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F~STATISTICS FOR THE 25. 7777۳ 

SUM OF SUMS = 0.8440126623444366D 02 

SUM OF SUMS**2 = 0.1782254895437147D 03 

SUM OF SQUARE OF SUMS = 0.6892901073945840D 03 

AMCNG GROUP SUM OF SQUARES = 0.2143851912627546D 01 

WITHIN GROUP SUM OF SQUARES = 0.17133258 84697689D 03 

TOTAL SUM OF SQUARES = 0.1734764403823964D 03 

DEGREES OF FREEDOM IN NUMERATOR - 0.1400000000000000р 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC - 0.1327251499840134D 01 


F-STATISTICS FOR THE 26 COMPONENT 

SUM OF SUMS = 0.6923579472492147D 02 

SUM OF SUMS**2 = 0.1481985058366125D 03 

SUM OF SQUARE OF SUMS = 0.4452980334 185913D 03 

AMCNG GROUP SUM OF SQUARES = 0.1257250153391605D 01 

WITHIN GROUP SUM OF SQUARES = 0.1437455255024266D 03 

TOTAL SUN OF SQUARES = 0.1450027756558182D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.92773892874 76910D 00 


F-STATISTICS FOR THE 27 COMPONENT 

SUM OF SUMS = 0.6366092822129744D 02 

SUM OF SUMS**2 = 0.1177906819912730D 03 

SUM OF SQUARE OF SUMS = 0.3935159041535181D 03 

AMCNG GROUP SUM OF SQUARES = 0.1233349853537058D 01 

WITHIN GROUP SUM OF SQUARES = 0.1138555229497378D 03 

TOTAL SUM OF SQUARES = 0.1150888728032749D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.1149027974258136D 01 
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F-STATISTICS FOR THE 28 COMPONENT 

SUM OF SUMS = 0.5387592494868525D 02 

SUM OF SUMS**2 = 0.1186133442790129D 03 

SUM OF SQUARE OF SUMS = 0.2792595957915495D 03 

AMCNG GROUP SUM OF SQUARES = 0.8575190985312517D 00 

WITHIN GROUP SUM OF SQUARES = 0.1158207483210974D 03 

TOTAL SUM OF SQUARES = 0.1166782674196286D 03 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D 04 
F-STATISTIC = 0.7853366268738306D 00 


F-STATISTICS FOR THE 29 COMPONENT 

SUM OF SUMS = 0.4536144151906700D 02 

SUM OF SUMS**2 = 0.6564305139058856D 02 

SUM OF SQUARE OF SUMS = 0.1919053830677748D 03 

AMONG GROUP SUM OF SQUARES = 0.5472802462192574D 00 

WITHIN GROUP SUM OF SQUARES = 0.6372399755991081D 02 

TOIAL SUM OF SQUARES = 0.6427127780613007D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0. 1485000000000000D 04 
F-STATISTIC = 0.9 109723144851782D 00 


F-STATISTICS FOR THE 30 COMPONENT 
SUM OF SUMS = 0.3732568029970994D 02 

SUM OF SUMS**2 = 0.5705765228634419D 02 

SUM OF SQUARE OF SUMS = 0.1284820507967475D 03 

AMCNG GROUP SUM OF SQUARES = 0.3560162347433721D 00 

WITHIN GROUP SUM OF SQUARES = 0.5577283177837672р 02 

TOTAL SUM OF SQUARES = 0.5612884801312009D 02 

DEGKEES OF FREEDOM IN NUMERATOR - 0.14000000000000005 02 
DEGREES OF FREEDOM IN DENOMINATOR - 0.1485000000000000D 04 
F-STATISTIC = 0.6770886363437515D 00 
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Ғ-5ТАТІ5ТІС5 ҒОК ТНЕ 31 COMPONENT 

SUM OF SUMS - 0.2973876286269883D 02 

SUM OF SUMS**2 - 0.42178336 10659068D 02 

SUM OF SQUARE OF SUMS - 0.7726205204918280D 02 

AMONG GROUP SUM OF SQUARES = 0. 18302 45094226048р 00 

WITHIN GROUP SUM OF SQUARES = 0.4140571558609885D 02 

TOTAL SUM OF SQUARES = 0.4158874009552145D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0-.14850000000000000 О4 
F-STATISTIC = 0.4688645251806333D 00 


F-STATISTICS FOR THE 32 COMPONENT 

SUM OF SUMS - 0.2142367689889270D 02 

SUM OF SUMS**2 = 0.1414820425416038D 02 

SUM OF SQUARE OF SUMS = 0.3565723402182890D 02 

AMONG GROUP SUM OF SQUARES = 0.5058971897285643D-01 

WITHIN GROUP SUM OF SQUARES = 0.1379 163191394209D 02 

TCIAL SUM OF SQUARES = 0. 13842221632914S55D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.14850000000000005 04 
F-STATISTIC = 0.38908548284%65455D 00 


F-STATISTICS FOR THE 33 COMPONENT 
SUM OF SUMS = 0.1732004830070321D 02 

SUM OF SUMS**2 - 0.111558209706648 1D 02 

SUM OF SQUARE OF SUMS - 0.2312899538600082D 02 

AMCNG GROUP SUM OF SQUARES = 0.3130057176754676D-01 

JITHIN GROUP SUM OF SQUARES = 0.1092453101680480D 02 

TOTAL SUM OF SQUARES = 0.1095583158857235D 02 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.1485000000000000D O4 
F-STATISTIC = 0.3039120267386334D 00 
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F-STATISTICS OF MEAN VARIANCE, SKEWNESS AND KURTOSIS 


F-STATISTICS FOR THE 1 COMPONENT 

SUM CF SUMS = 0.1111617439422250D 05 

SUM OF SUMS*¥*2 = 0.4319425777749392D 07 

SUM CF SQUARE OF SUMS = 0.18612607339607520 08 

AMONG GROUP SUM OF SQUARES = 0.3143833877400904D 06 

WITHIN GROUP SUM OF SQUARES = 0.3755407373518861D 07 

TOTAL SUM OF SQUARES = 0.4069790761258952D 07 

DEGREES OF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = O.480000 00000000000 03 
F-STATISTIC = 0.2870223636519026D 01 


Pe Lat otiGs FOR THE 2 COMPONENT 

SUM OF SUMS =-0. 3062348822555926D 07 

SUM CF SUMS**2 = 0.7751235739069674D 12 

SUM OF SQUARE OF SUMS = 0.2376364181020636D 13 

AMONG GROUP SUM OF SQUARES = 0.5306562101878761D 11 

WITHIN GROUP SUM OF SQUARES = 0.70311253811846348D 12 

TOTAL SUM OF SQUARES = 0.7561781591372509D 12 

DEGREES OF FREEDOM IN NUMERATOR = 0.14000000000009000 02 
DEGREES CF FREEDOM IN DENOMINATOR - 0.48000000000000060D 03 
F-STATISTIC = 0.2587626620217669D 01 


Р-5ТАТЇ511С5 FOR THE 3 COMPONENT 

SUM OF SUMS = 0.10876927598750710 08. 

SUM OF SUMS**2 = 0.1368235571001090D 14 

SUM OF SQUARE OF SUMS = 0.3521258942662316D 14 

AMCNG GROUP SUM OF SQUARES = 0.8280430048704694D 12 

WITHIN G&OUP SUM OF SQUARES = 0.1261530754556778р 14 

TOTAL SUM OF SQUAKES = 0. 1344335055043825D 14 

DEGREES CF FREEDOM IN NUMERATOR = 0.1400000000000000D 02 
DEGREES OF FREEDOM IN DENOMINATOR = 0.4800000000000000D 03 
F-STATISTIC - 0.22504442140596077D 01 
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